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NAVAL ARCHITECTURE 


Stability of free-floating ship 


Part I 


Maciej Pawtowski 
Gdansk University of Technology 


ABSTRACT 


Problem of calculation of righting arms of the free-floating ship, i.e. longitudinally balan- 
ced at any heel angle, was formulated. In such a case of particular interest for a ship in 
the damage condition, the righting arms are ambiguous as they depend on a way the 
heeling moment acts. Two cases were considered : when the heeling moment is parallel to 
the ship plane of symmetry, and the case when it performs the least work, i.e. when the 
moment is parallel to the main axis of ship waterplane. It was demonstrated that angular 
translations (heel and trim) are then the Euler angles associated with a relevant reference 


axis. Some cases of the incorrect defining and using of those angles in today design practice were indica- 
ted. The most important features of the curve of righting arms of free-floating ship were demonstrated. 


Keywords : stability, free floating ship 


INTRODUCTION 


Rules of the classification societies require the intact ship 
stability to be investigated for the ship floating at even keel. 
The rules, however, do not clearly state how to calculate the 
damaged ship stability, which often leads to large discrepances 
in obtained results. 

For the undamaged ship it is practically meaningless whe- 
ther the stability calculations are performed for the ship ha- 
ving a fixed trim, constant in function of the heeling angle, or 
for the free-floating ship which changes its trim depending 
on its longitudinal equilibrium state. This is due to a small 
asymmetry of the ship relative to its midship section plane. 
However for the damaged ship the problem is important as it 
significantly influences the course of the righting arm curve 
for the heeeling angles greater than the entrance angle of the 
deck into water (Fig.1). The righting arm is understood here 
as the distance between the action line of buoyancy force and 
that of gravity force, occurring in still water, at a given he- 
eling angle. 

The influence of calculation assumptions is especially im- 
portant in the case of flooding compartments far off the mid- 
ship, which is understable due to the then occurring high asym- 
metry and a small entrance angle of the deck into water. More- 
over the influence very strongly increases along with the ratio 
L/B decreasing. Hence it is greater for catamarans and SWATH 
units. The differences between righting arms may even reach 
a few hundred percent. For this reason the regulations should 
clearly define a way of carrying out calculations : whether at 
a fixed or free trim. It should be remembered that the final aim 
of stability calculations is to determine an expected final state 
of a considered ship under action of disturbing moments — and 
as a result — to correctly assess stability safety of the ship. 
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Fig. 1. Righting arm curves for a free-floating damaged ship 
and that having a fixed trim [1] 


It is obvious that the routine stability calculations should 
be performed for a free-floating ship. However in such com- 
monly accepted case the problem of correctly understood ship 
heel angle arises because then it is an ambiguous notion, which 
is manifested in existing various definitions of that angle. Most 
often the heel angle is assumed to be the slope angle of water- 
-level trace line on the frame planes, further marked @, or the 
slope angle of the baseplane (BP) relative to water-level , fur- 
ther marked a. The first of the angles is loosely associated with 
correctly defined ship heel angle, whereas the second is the 
heel angle of the ship having minimum righting arms. 


HEEL ANGLE OF FREE-FLOATING SHIP 
The problem of correct defining the heel angle of the free- 


-floating ship has been recently solved [2+4]. Namely, the ship 
heel angle, further marked þ, is understood as the angle rota- 
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tion of the ship plane of symmetry (PS) around the intersection 
line of water surface and PS, i.e. the angle of rotation around 
the water-level trace line on PS, in other words — this is the 
inclination angle of PS from the vertical. At the same time the 
angle is equal to the slope angle of the axis y relative to water 
surface. The definition stems from the assumption that the ship 
heels under action of the heeling moment of fixed direction in 
space, and parallel to PS. If the ship floats freely it heels in such 
a way as to be logitudinally balanced all the time. It means that 
the direction of the righting moment is fixed in space - and the 
same applies to the heeling moment. The conclusion immedia- 
tely follows that the curve of buoyancy centres is then exactly 
flat and situated on the rotation plane perpendicular to the he- 
eling moment vector, and containing the ship centre of gravity. 

A necessity to define the heel angle of free-floating ship is 
usually not felt — many surveyors and designers are just surpri- 
sed that any problem of this kind exists at all. Hence various 
definitions of the angle in question have been still assumed, 
which obviously results in ambiguity of calculation and makes 
it not possible to compare different computer softwares. 


HISTORICAL OUTLINE 


Why a body floats in a liquid has been known already in 
antiquity since the times of Archimedes. However in which way 
to assess and investigate stability of floating bodies has been 
recognized only after discovery of the Newtonian laws. In 1746 
Bouguer introduced the notion of metacentrum and metacen- 
tric height considered as a measure of initial stability. In 1749 
Euler introduced a formula for metacentric radius, and a theo- 
rem for equi-volume waterplanes. In 1796 Atwood published 
a method for calculation of the righting arm at a given heel 
angle [5]. Nonetheless for over a hundred years only the initial 
metacentric height GM was used to assess ship stability. Only 
the ship stability accidents at the end of 19' century revealed 
importance of the freeboard and necessity of applying the cur- 
ve of righting arms in assessing stability of ships. 

The metacentric height, which is not an unimportant index 
of stability, does not make it possible to directly assess either 
a range of the curve of righting arms nor on a value of the 
maximum righting arm. Here it is worth mentioning the wide- 
ly described case of sinking the HMS Captain in 1870, whose 
metacetric height GM = 0.79 m [6]. The ship capsized during 
a storm in the Bay of Biscay, whereas the accompanying bat- 
tleship Monarch of a similar size and characteristics survived 
unharmed despite having its metacentric height GM = 0.73 m, 
i.e. smaller than that of the first ship. The fact was very surpri- 
sing for the then naval architects. It is very easy to explain the 
accident if one observes that the freeboards of the two ships 
much differred to each other : the Captain had the freeboard 
F = 1.98 m and the Monarch - F = 4.27 m. As a result, despite 
the smaller metacentric height of the Monarch, its curve of 
righting arms was of much better parameters than that of the 
Captain, whose GZmax = 0.55 m instead of 0.25 m, Omax = 40 
deg instead of 19 deg, and the stability range angle 0, = 70 deg 
instead of 54 deg. 

The Captain’ accident has finally proved that the meta- 
centric height is an insufficient measure of safety against cap- 
sizing and it has made it necessary to examine ship’s stability 
also at large heel angles. As a result, at the end of the 19° 
century the curves of righting arms, called the Reed 5 curves in 
memory of their propagator, began to be used for the ship sta- 
bility assessment. The first stability criteria, given by Rahola 
[7], appeared as late as in 1939. The well-documented recom- 
mendations dealing with a minimum size of the curve of righ- 
ting arms have been elaborated on the basis of the analysis of 
the curves of righting arms of both for capsized ships and sta- 
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ble ones. At the end of the 1960s those criteria were adopted 
by IMCO (Intergovernmental Maritime Consultative Organi- 
zation, presently IMO (International Maritime Organization), 
and they have been valid until now [8]. 

Though the curve of righting arms has been applied to as- 
sess stability of intact ships, stability of damaged ships has been 
further controlled by means of the metacentric height and fre- 
eboard. In the SOLAS conventions including the last one of 
1974 the residual freeboard of as low as only 3 inches and the 
metacentric height of 2 inches have been assumed permissible. 
With such parameters the curves of righting arms usually show 
marginal values. A change took place as late as in 1990 when 
criteria for the curve of righting arms were introduced for da- 
maged ships, in the form of the SOLAS 90 criteria [9]. How- 
ever it is worth remembering that the criteria have not repre- 
sented any important progress as they resulted from an admi- 
nistrative decision. Therefore they have only an alleged, but 
not real, relation to actual safety of a ship in damage condition. 
A breakthrough in that regard has happened during the last six 
years [10,11]. 


FORMULATION OF THE PROBLEM 


Almost all of the commonly known calculation methods of 
the curve of righting arms concern the ship floating on even 
keel. This means indirectly that the buoyancy centre is assu- 
med not to translate longitudinally during ship’s heeling. It has 
been no necessity to consider earlier a different situation as the 
calculations dealt only with intact ships for which the assump- 
tion has been well satisfied. However the fact cannot be ne- 
glected any longer in the situations when the buoyancy centre 
translates longitudinally due to an asymmetrical distribution 
of buoyancy relative to the plane of rotation, as in the case of 
semi-submersible platforms arbitrarily orientated relative to 
wind direction, ships of low L/B ratio, or ships in damage con- 
ditions, and then the calculations should be carried out for 
a free-floating object. Determination of the curve of righting 
arms becomes in such cases ambiguous and the problem must 
be defined. Especially the way in which the righting moment 
acts should be defined. 

It should be said that angular motions of a free-floating ob- 
ject are not considered in the basic ship theory as this is a 3-D 
motion. And, such motion is spatial and requiring good spatial 
imagination. For this reason as well as for making calculations 
easier the vectorial calculus is applied in this paper. 


Calculations of the righting arms curve for free-floating 
ship is carried out under the following assumptions : 


a) a pure heeling moment is statically exerted on the ship. 
It means that ship inclinations are equi-volumetric and the 
horizontal position of the ship gravity centre remains con- 
stant (the moment cannot induce any translational motion 
hence change location of ship gravity centre on the sea le- 
vel) 

the heeling moment vector is strictly horizontal. If this is 
not the case a vertical component of the moment, able to 
rotate the ship around its vertical axis, will exist 

c) the heeling moment direction is fixed in space and paral- 
lel to PS. As the moment is simultaneously parallel to the 
water-plane hence it is also parallel to the trace line of wa- 
ter on PS. Therefore the intersection edge of PS and the 
water-plane is also fixed in space. The edge determines 
orientation of the ship at the sea level as well it defines the 
direction of the heeling moment ”) 

the ship is in static equilibrium, i.e. the sum of forces and 
moments acting on it equals zero. Hence ship’s weight is 
equal to its buoyancy, and the statically applied heeling mo- 


b 


— 


d 


— 


e) 


ment is balanced by the righting moment of the same direc- 
tion and opposite sense 

the righting moment is formed by the couple of forces : 
i.e. the gravity force applied in the ship’s centre of gravity 
and the buoyancy force applied in the ship’s centre of bu- 
oyancy — the forces are equal and of opposite sense to each 
other. The moment’s vector is horizontally directed, per- 
pendicular to the vertical plane determined by the gravity 
force and buoyancy force. 


From those assumptions some consequences follow : 


As in the state of equilibrium directions of the moments are 
the same the centre of buoyancy must be situated on the 
plane perpendicular to the direction of action of the heeling 
moment (i.e. on the rotation plane), and on which the cen- 
tre of gravity is located. As the direction of the heeling 
moment action is, under the assumption, fixed in space, 
hence the rotation plane is also fixed in space. The versor 
perpendicular to the rotation plane, further marked e, is 
hence constant and stands for the rotation axis. 


For ship heels at a fixed trim the centre of buoyancy need 
not to be located on the rotation plane, therefore the mo- 
ment acting on the ship has not a constant direction at the 
horizontal plane. 


The rotation plane is vertical, motionless in space, perpen- 
dicular to the direction of the moment action (trace of wa- 
ter on PS), and passing through the ship gravity centre G 
motionless at the horizontal plane. 


It is the rotation plane about which the buoyancy centre of 
free-floating ship moves. Hence the curve of buoyancy cen- 
tres is strictly flat in space. For ship heels at a fixed trim the 
curve of buoyancy centres is the projection of the spatial 
curve onto the rotation plane. 


The trace of water on PS determines the direction of the he- 
eling moment. The ship heel angle ọ is the rotation of PS 
around the trace of water on PS, i.e. the deflection angle of 
PS from the vertical. As Oy-axis is perpendicular to PS the 
angle ¢ is simultaneously the slope angle of that axis against 
the water-level. The angle ọ can be interpreted as the rota- 
tion angle of Oy-axis around the trace of water on PS, in the 
rotation plane. 


In order the ship to be balanced at any heel angle it has 
also to rotate around the versor normal to PS, i.e. the Oy- 
-axis, in such a way as to bring the centre of buoyancy on 
the rotation plane (Fig. 2). The ship rotation angle in PS, 
further marked 9, is called the trim angle — this is the defi- 
nition commonly used in ship hydrostatics. A change of 
the trim angle O does not influence the heel angle ọ. The 
two angles are the Euler angles associated with Oy-axis. 
The third angle associated with the rotation of Oy-axis pro- 
jection onto the horizontal plane (i.e. yaw) is not present as 
the rotation plane is fixed in space. Under the above given 
assumption b) any yaw-inducing moments are not conside- 
red in ship hydrostatics. 


The righting arm GZ is the arm of the couple of forces for- 
ming the righting moment, i.e. the distance, measured in 
the rotation plane, between the action line of gravity force 
and that of buoyancy force. 


As the righting moment is all the time parallel to the trace 
of water on PS, the righting moment work is the integral of 
the moment, respective to the heel angle 6. Simultaneously 
this is the least work which is to be performed in order to 
heel the ship up to a given angle @ . In other words, for 


a ship of a fixed trim or that not fully balanced the work of 
righting moment is greater. 


Fig. 2. Rotation plane of free-floating ship which trims in PS 


It can be observed that the projection of Oy-axis onto the 


horizontal plane is perpendicular to the trace of water on PS. 
Hence the presented model of the heeling moment action stric- 
tly corresponds to the heeling moment due to a shift of cargo in 
ship’s transverse plane. This deals also with the heeling mo- 
ment of ro-ro ships in the damage condition, resulting from the 
accummulation of water on the car deck when a symmetrical 
midship compartment has been flooded. For the same reason 
the curve of righting arms measured by means of the Di Belli 
method?) is strictly consistent with the above given model of 
inclinations. Finally, the heeling moment of fixed direction in 
space, parallel to the trace of water on PS, is an accepted ideali- 
zation of the wind-induced moment. 


BASIC RELATIONSHIPS 


A right hand side coordinate frame Oxyz, shown in Fig.3, 


fixed to the ship, is assumed. The point O is identical with the 
point K, Oy-axis points port, and Oz-axis — upwards. 


Fig. 3. The right hand side coordinate frame Oxyz 


An arbitrary attitude of the waterplane 
can be described by the equation : 


z= To + xtg0 + ytgo (1) 


in which three independent parameters appear. These are : 


+ the angle 0 called the trim angle which is the slope angle 


of the trace of water on PS relative to Ox- axis 
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+ the angle @ which is the slope angle of the trace of water on 
the midship plane relative to Oy-axis, and 
+ T,- draught of Oz-axis. 


The above mentioned angles are positive when a positive 
increment of z corresponds to a positive increment of x (or y) , 
as shown in Fig.3. Therefore the trim angle O > 0 is positive 
when the ship is trimmed by bow, and the angle @ > 0 is posi- 
tive when it is heeled a-port. Both the angles are easy to be 
measured as : 


tg0 = ATpr/Lpp and tg@ = AT, p/B 


where : 


ATpr — difference of bow and stern draughts 
measured at the respective perpendiculars 
ATıp — difference of port and starboard draughts 


measured midships. 


The angles @ and @ unambiguously describe position of the 
ship against water-level (or water-level against the ship) and 
are called the analytical angles. From the analytic geometry it 
results that the vector R normal to the waterplane given by (4) 
has the components : 


R= (tg0 > tgọ > -1) (2) 


which points downwards and its absolute value equals : 


R =4/1+tg’0+tg’ọ (3) 


Hence the versor n normal to the waterplane 
and pointing upwards is equal to : 


=—-R/R (4) 


The angle & contained between the water-level and the base 
plane BP (or initial waterplane) is one of the heel angles taken 
for calculation of the curve of righting arms. The angle betwe- 
en the planes is the same as that between the versors n and k 
normal to them. As a result, cosa = k - n. Therefore : 


cosa =1/R=1/ /1+tg’0+ tg’ 


Hence : 


tga /tg°0+ tg’ (5) 


The sign of the angle is the same as that of the angle @. 
Taking into account that 1/R = cosa the formula (4) yields the 
following components of the versor n : 


n = (-tgOcosa , -tg@cosa , cosa) (6) 


The heel angle is equal to the inclination angle of PS from 
the vertical. The angle is the same as that of Oy-axis relative to 
water-level. Hence cos (90° + 6) = j > n = ny. And, sing = 
= tg@cosa, or simpler : 


tgb = cosOtg@ (7) 


The rotation plane rotates around the rotation axis paral- 
lel to the direction of the heeling moment action, defined by 
the versor e. The rotation angle is equal to the angle ọ . The 
heeling moment direction is parallel to the water-level trace 
line on PS, whose versor e; = (cos0,0,sin@), see Fig.3. Hence 
e= e]. 

Righting arm 


In the course of heeling the ship its displacement remains 


constant and its buoyancy centre shifts in the rotation plane 
perpendicular to the rotation axis. Therefore the following is 


yielded : 
e-r=0 (8) 
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for any heel angle , where : 


r = (XB — XG , YB — YG; ZB — ZG) — buoyancy centre radius 
vector relative to the ship gravity centre r = GB. 


The righting moment is given by the formula : 
M=r-nA (9) 
where : 
A = gV - ship buoyancy. 
The vector is parallel to the rotation axis e, hence : 
M=e-(r-nj)A 
Therefore the righting arm GZ = M/A 
is given by the formula : 
GZ=e: (r:n) 


which is a function of the heel angle defined by (7). It is 
worth mentioning that the angle  < a, which can be observed 
directly from the formula cosa = cosOcos obtained by divi- 
ding sind by tgb. From (7) it also results that ġ < @. Therefore 
the actual rotation angle ¢ is never greater either than the angle 
a or the angle @, which have been taken as the heel angles of 
free-floating ship, so far. 


(10) 


Metacentric radii 


The buoyancy centre of the free-floating ship moves along 
a curve in the rotation plane which rotates as a disc around the 
rotation axis motionless in space. As the action lines of the 
buoyancy force are always vertical they are perpendicular to 
the actual waterplane. When changing the ship heel by dọ the 
buoyancy force lines rotate also by the angle dọ in the rotation 
plane, and their respective waterplanes — by the angle do, aro- 
und the instantaneous axis of waterplane rotation, called the 
ship floatation axis f. The relationship between the differen- 
tials is given by the following formula, [4 ] : 


da cosx = do (11) 
where : 
x, -the angle defining situation of the floatation 


axis relative to the rotation axis e = e1, 
being the PS trace line on waterplane. 


The equation (11) reflects the fact that small angles have 
vectorial features. Hence the angle dọ is nothing else but the 
projection of the waterplane rotation angle da, onto the rota- 
tion axis e. In general case the angle do, is not equal to the 
change of the slope angle da of waterplane relative to BP; the 
equality occurs only when the floatation axis f is parallel to the 
intersection edge of BP and water-level (see Fig.3). 

The metacentric radius is understood as the curvature ra- 
dius of the curve of buoyancy centres in the rotation plane; 
it depends on the angle ©. Translation of the buoyancy centre 
along the arc of the curve of buoyancy centres amounts to ds = 
= redo, = BMdọ. On accounting for the formula (11), the meta- 
centric radius can be determined by the formula : 


Tp 
COSY, 
where : 
n =(1/VW I} +D? 
f — instantaneous ship’s floatation axis passing 


through the floatation centre F (the waterplane 
centre of gravity) 
V — ship’s hull volumetric displacement 


I¢ and Df — transverse and deviation (cross-product) iner- 
tia moment of the instantaneous waterplane, 
respectively, associated with the floatation axis 
and floatation centre F. 


The quantity rg is a proportional coefficient 
between the buoyancy centre translation ds and the angle da. 


The centre of buoyancy translates in the rotation plane, in 
parallel to the instantaneous waterplane (water-level). There- 
fore the buoyancy centre translation vector dr = (n - e)ds. 


Floatation axis 


As it can be observed the metacentric radius of the free- 
-floating ship , at a given heel angle, depends on location and 
orientation of the instantaneous floatation axis f. If the water- 
plane is rotated by the angle do the transverse component of 
the buoyancy centre translation, BC, relative to the floatation 
axis (Fig. 4) is proportional to If, and the longitudinal compo- 
nent of this translation, AB,— to Ds, which results from the 
Euler’s theorem of equi-volume waterplanes; this is the reason 
that the expression for the resultant translation AC, equal to 
ds = rpdoy , appears in (12). 


Fig. 4. Waterplane top view 


It is required the resultant translation to be normal to the 
direction of heeling moment action (the rotation axis e). In 
order to obtain this the angle C in Fig.4 must be equal to x, 
which results from the properties of the angles having corres- 
pondingly perpendicular arms. Therefore the floatation axis 
slope angle relative to the rotation axis has to satisfy the follo- 


ion: 

wing equatio tey = Dr ie (13) 

The angle x has the same sign as that of the waterplane 
deviation moment (in Fig.4 - positive). It should be remembe- 
red that the moments Dr and I¢ also depend on x. By applying 
the relationships resulting from the Mohr circle [4,12], known 
from the theory of strength of materials, they can be represen- 
ted as follows: Dre= D', and If = I + a'. The equation (13) obta- 
ins then the form : D' — (I + a')tgy = 0, where the mark ,,'” 
stands for the quantities associated with the rotated coordinate 
system (&', n') whose axis €' coincides with the floatation axis 
f, and its origin is located in the floatation centre F; (the system 
is not shown in Fig. 4). The quantities marked ,,'” are expres- 
sed by means of the quantities taken from the system (€, n) : 


D' = Deos2y, + asin2x = rsin(2y + 2%) 
a' = acos2x — Dsin2x = rcos(2y + 2%) 
where : a = '(Ig¢ — Inn)! = “(lee + Inn) and 
D = len, lee, Inn - waterplane deviation moment, trans- 
verse and longitudinal inertia moment, 
respectively, in the system én (Fig.4) whose origin coinci- 


des with the floatation centre F, and the axis € is parallel to 
the rotation axis e = e, (a trace of PS on the waterplane). 


The quantities D' and a' determine the parametric equation 
of the Mohr circle shown in Fig.5. When applying the introdu- 
ced notation the equation (13) gets the following form : 


rsin(2y + 2x) — [I + rcos(2y + 2x)]tgy =0 (14) 
where : 
r= va’+D’ , 2y, =arctg(D/a) 


2y = 2yo if a > 0, in the opposite case : 2y = 2y, + 180° 


The equation (13) having the quantities marked ,,'”’ is easier 
to be solved, and the equation (14) more simple for the geome- 
trical interpretation shown in Fig.5; where a and Yo are assu- 
med negative. 


Inertia 
moments I 


Deviation moments D 


Fig. 5. Mohr circle and geometric characteristics of waterplane 


When cos2x and sin2x appearing in (13) is expressed by 
tex the equation (13) can be reduced to the simple 1% order 


tion : 
equation D = (1—a)tgy (15) 
Therefore : tg% = D/Inn 
where : |x| < arcsin(r/I) — see Fig.5. 


As the angle ¥ is known, the quantities De= D' and Ip=I+ a’, 
necessary to express the radius on the basis of (12), as well as 
direction of the floatation axis, are defined. The floatation axis 
versor is given by the formula : 

f = ecos% + (n- e)siny (16) 

The equation (13) results from the asssumption that : e - dr = 
= 0, i.e. that the buoyancy centre translation in the ship-fixed 
coordinate system is normal to the rotation axis. It would be 
this way if the rotation axis were fixed to ship; however the 
axis is fixed in space but this does not mean it is fixed to the 
ship. It can be observed that when changing the heel also, in 
general, the trim has to be changed to balance the ship, which 
results in changing the orientation of the axis e; = (cos0 , 0, 
sin@) relative to ship, determined by a trace of water on PS. 

After differentiation of the equation (8) the following is 
obtained : e: dr =-— de -r, i.e. that the buoyancy centre transla- 
tion in the hull-fixed coordinate frame is not strictly normal to 
the rotation axis. It should be intuitively obvious : the buoyan- 
cy centre translation in the ship-fixed frame has to be oblique 
to it as the buoyancy centre is to be permanently located in the 
rotation plane which changes its orientation relative to ship 
during inclinations. When this is accounted for, the following 
relationship for the angle x between the floatation axis and 
rotation axis is obtained [4] : 


D 
tgx = = (17) 
I,-BZV V(BM, -BZ) 
where : 
BZ =-r- n —is the height of the gravity centre over the buo- 
yancy centre (Fig.2) 


BM, — is the longitudinal metacentric radius repre- 


sented by the bracketed term in (17). 


As the term BZV is negligibly small in comparison with the 
longitudinal inertia moment of waterplane, Inn, the equation 
(17) practically yields the same solution as the equation (15). 
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Hence it can be observed that to determine the longitudinal 
radius of free-floating ship at a given heel angle it is necessary 
to know three geometrical characteristics of waterplane, viz. 
its deviation moment and transverse and longitudinal inertia 
moments in the En coordinate system associated with the rota- 
tion axis e = e, (a trace of PS on the waterplane) and the buoy- 
ancy centre. In calculating these characteristics it should be 
accounted for that the water trace lines on frame planes are not 
normal to the water trace line on PS. Denoting this angle by 
90° + P the following is obtained : cos(90° + ) = e1 > e2, 
where : e2 = (0 , cos@ , sin@) is unit versor of the trace of water 
on frame planes (Fig.3). 

Hence : 
sinB = —sinOsing (18) 

The way of calculation of the geometrical characteristics 
of the waterplane arbitrarily heeled relative to ship is discussed 
in [4] and [13]. 

When the floatation axis f is known it is easy to find the 
analytic angles @ and O describing ship’s location relative to 
water at a new heel angle. Namely, by changing the slope of the 
waterplanes relative to each other, by the angle Ac, , the rota- 
tion of the versor n around the floatation axis by the angle Aq), 
is induced. Hence the new versor n4 is given by the formula : 


n; = ncosAq, + (f - n)sinAq, (19) 


As the new versor n is known, the new analytical angles 
corresponding to the versor can be easily found by using the 
formula (6). Namely, tg0 = —nx/nz, and tg@ = —ny/nz. Because 
the translation dr and the new waterplane slope angles are 
known, it is very fast to find a balanced location of buoyancy 
centre at the new heel angle enlarged by the angle Ao . 

The versor n can be directly expressed by the angles O and 
Q serving as degrees of freedom. To this end, in the formula (6) 
the expression cosa = cosOcos@ as well as that for tg given 
by (7) should be accounted for. The following is immediately 
obtained : 


n = (-sinOcosd , -sinh , cosOcos) (20) 


Knowing the new versor n, one has as before : tg9 = —n,/nz, 
and sino = —ny. It is not necessary to find tg@ =—ny/nz, as in the 
first approximation the equation of new waterplane (in the ship- 
-fixed frame) can be represented by the formula : 

ny(X — Xp) + ny(y — yg) + n,(z - zp) = 0 
where : 
Xp, yr, Zp - the coordinates of the previous floatation centre. 


(21) 


The equation (21) is more convenient than (1) , as tg@ and 
To tend to infinitely large values along with the heel angle in- 
creasing. The equation (1) is necessary to start calculations. 
The waterplane given by (21) is neither equi-volume nor ba- 
lanced one. Correct values of draught and trim can be found by 
means of the method of successive approximations so as to 
maintain ship displacement constant and equal to an assumed 
value, and to make the equation (8) satisfied. The calculations 
become significantly shorter in case of making use of the pro- 
perties of equi-volume waterplanes of free-floating ship. 
Within a finite interval of the heel angle Ad , such waterplanes 
roll on the surface of a cone whose parameters can be determi- 
ned in advance [4]. The rolling waterplanes adhere to the cone 
along an instantaneous floatation axis. 


Mechanism of equi-volume heels 


An infinitesimal rotation of waterplane around the floata- 
tion axis f can be considered as resulting from two rotations : 
ship’s rotation by the angle dọ around the rotation axis 
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e = e; (trace of water on PS), and ship’s rotation by the angle 
dð around the normal to PS. The directed angle jd@ is inclined 
by the angle @ relative to water-level. As small rotations have 
vectorial properties the directed angle fda, is the resultant of 
two components : parallel and normal to the rotation axis e : 


fda, = (do , — dOcos) (22) 


The components are defined in the En coordinate frame 
(see Fig.4). The normal component to the rotation axis, equal 
to — dOcos@ is further denoted daz. The directed angle fda is 
inclined by the angle x to the rotation axis (€-axis). It can be 
seen in Fig.4 that the positive normal component da corres- 
ponds to the positive angle x, and the trim change is negative 
(by stern), therefore the normal component of the opposite sign 
must be applied. The projection of dœ; onto the rotation axis 
yields from the relationship (11). Taking into account the rela- 
tionships inherent to rectangular triangles, one can determine 
the normal component day in three different ways : 


doz = da siny = dotgx = — dOcoso (23) 


From this formula it results that : 


> the more inclined the floatation axis from the rotation axis, 
the greater changes of ship trim during heeling, which is 
consistent with intuition 


> when x = 0, i.e. f =e, ship trim does not change due to the 
waterplane rotation itself as in the case of fixed-trim ship 


> for p = 90° (PS is then horizontal), % = 0, i.e. the rotation 
axis (which does not then exist in the sense of trace of wa- 
ter on PS) determines the floatation axis f. 


The rotation angle of PS around the normal , equal to jdé , 
has also the vertical component : — d@sing. It makes the angle 
x, changing as a result of a change of orientation of the rotation 
axis (relative to ship hull), which occurs in the course of trim- 
ming. The change of the angle x, which results from the chan- 
ge of orientation of the rotation axis only, is further denoted by 
dXo. Hence the angle dy, =— dOsind. Accounting for the equa- 
tion (23), one obtains : 


dxo = da sinxtgh = dotgxtgo (24) 
If for a new waterplane the floatation angle y changes by 
dx the new floatation axis rotates against the previous one by 
the angle dy¢ = d% — dX, equal to the difference of both the 
changes. When the angle dx¥¢> 0 is positive then the new floa- 
tation axis f shifts towards the heel , i.e. it departs from the 
rotation axis. 
The equi-volume waterplanes roll on a cone whose axis is 
inclined relative to them by the angle € determined by the fol- 
lowing formula : 


tge = dy¢/doy 


To find the formula is very simple : it is enough to observe 
that tge is the ratio of the radius and the generatrix of the cone. 
When the angle dxf > 0 is positive the cone is located above 
the waterplanes, otherwise — below them. The vertex of the 
cone is located in the distance /, of the generatrix from the 
floatation centre F, given by 1 = — dn'p/dx¢, where dn'r is the 
translation of the floatation centre perpendicular to the floata- 
tion axis (when | > 0, the vertex is located fore). Taking into 
account that dn'p = rpdoone obtains : 


= — rpdo,/dy- = — rp/tge 


where : 


rp = dl¢/dV — differential metacentric radius 
(curvature radius of the curve of floatation centres). 


From the formula it results that the cone base radius at the 
level of the floatation centre equals the differential metacentric 
radius. 

The normal waterplane rotation component dc represents 
the ship trim angle measured in the vertical plane passing thro- 
ugh the trace of water on PS. If ship heel is enlarged by do, the 
buoyancy centre translation perpendicular to the rotation angle 
is proportional to Dd, where : D - waterplane deviation mo- 
ment in the én coordinate frame (Fig.4) . 

The translation must be compensated by the trim Ippdap. 
Equating them to each other one gets da = (D/Inn)do. Accoun- 
ting for that tg% = do/d, one obtains the formula (15). A more 
exact solution can be obtained by using the metacentric formula 
for doy = (D/VGM_z )do, where : GM, = BM; — BZ is the longi- 
tudinal metacentric height. As tg% = do2/dd, the above gives 
the formula (17) provided before without derivation. The for- 
mula accounts for the rotation plane slope change defined in 
the hull-fixed coordinate frame, which results from the ship 
trim change. 

(to be continued) 


NOMENCLATURE 

a - height of gravity centre over buoyancy centre in 
upright position of ship 

B - buoyancy centre 

BM - transverse metacentric radius 

BM; - longitudinal metacentric radius 

BP - baseplane 

BZ - height of gravity centre over buoyancy centre 

e - direction of rotation axis (versor normal to rotation 
plane) 

ej, & - versors of trace of water on PS and on midship plane, 
respectively 

f - versor of floatation axis 

F - freeboard 

g - gravity acceleration 

G - ship gravity centre 

GM - metacentric height 

GZ - righting arm 

i, j, k - versors of the ship—fixed coordinate frame whose 


origin is in the point K (intersection point of the plane 
of symmetry, PS, midship plane and base plane BP) 


l, la, - righting arm and dynamic arm, respectively 

L; B, T - length, breadth and mean draught of ship, 
respectively 

n, - upward pointing versor nomal to waterplane 

PS - ship plane of symmetry 


SOLAS - Safety of Life at Sea 
SWATH - Small Waterplane Hull 


To - draught of z-axis 

WEGEMT - European Association of Universities in Marine 
Technology 

V - volumetric displacement of ship 

(0A - angle between initial waterplane and water-level 

B - angle between water trace lines on PS and midship 
plane 

A - ship buoyancy (weight of the displaced water) 

Al - correction of righting arm obtained by means of 


cross-curves of stability, accounting for oblique 

translation of gravity centre relative to rotation plane 

at changing the height of ship gravity centre over BP 

rotation angle of plane of rotation 

- slope angle of trace of water on PS relative to x-axis 
of ship 

- water density 

- slope angle of trace of water on the stations relative 

to y-axis of ship 

- angle of PS inclination from the vertical 

angle of vanishing stability 

- angle between floatation axis and rotation axis 

- angle between water-level trace line and PS trace line 
on initial waterplane 


os 
1 


SAFS sv 


D) Fix in space” does not mean : fixed relative to the ship hull coor- 
dinate system. 

In case of the objects arbitrarily situated relative to the heeling 
moment vector or wind direction (e.g. semi-submersible floating 
units) the comments concerning PS should be applied to the refe- 
rence plane initially situated perpendicularly to the wind direc- 
tion and rotating together with the object in question. For the ship 
unsymmetricaly flooded this is the plane parallel to the main in- 
ertia axis of the waterplane in ship’s upright position. 

In the Di Belli method, a heel angle of ship model, induced by 
shifting a weight along an arm perpendicular to PS, is measured. 
The heel angle is the inclination angle of the arm against water- 
-level. 
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ABSTRACT 


Combustion piston engine is one of the devices in which fast changeable processes occur in operational 

conditions. In this paper are presented basic problems associated with research on fast changeable proces- 

ses occurring in diesel engines , exemplified by the processes of indicated pressure and fuel pressure injec- 

ted to engines cylinder. Dynamical characteristics of the investigated processes were analyzed and pro- 

blems of synchronous averaging of pseudo-periodical signals were considered in order to limit high frequ- 

ency noise content in useful signal. Some limitations of elimination effectiveness of high frequency noise 
from tested signals have been revealed. 


Keywords : piston combustion engine, fast changeable processes, stochastic processes, high frequency noise 


INTRODUCTION 


Investigations of fast changeable processes, related to the 
dynamical features characteristics typical for operational con- 
ditions of devices, cause many methodological problems. This 
is especially important because of a high level of ignorance on 
disturbances of the processes associated with measurement 
systems as well as occurrence — in real objects — of factors not 
directly accounted for in testing program. The piston combu- 
stion engine is a typical object in description of which such 
problems appear [3 + 6]. 

The processes occurring in the combustion engine can be 
classified from the point of view of their dynamical features, 
and a role they play in the engine’s operation, namely [3] : 


> processes associated with particular working cycles of the 
engine as well as the fast changeable ones in comparison 
with those characteristic for the typical use of the engine 


> processes associated with the typical use of the engine 


> processes associated with service (tribological) wear of the 
engine as well as slowly changeable ones in comparison 
with the processes characteristic for the typical use of the 
engine. 


The first group of the processes contain all the processes 
connected with engine’s cycle. Processes not fully connected 
with the cyclic work of the engine, however generated by it, 
e.g. vibration of its elements, can be also considered as belon- 
ging to that group. Frequencies of the processes, conventio- 
nally considered as characteristic ’ for them, should be at le- 


ast one-order higher than those of engine working cycles as 
they occur in intervals corresponding to fragments of single 
cycles [3]. 

The processes characteristic for the typical use of the engi- 
ne are associated with operational forcing factors affecting the 
engine. The main service forcing factors are the following : 
control of the engine by its operator, and its loading from the 
side of power consumer. 

The power consumer reacts by rate of its task realization, 
and the engine — by changing its rotational speed. On the basis 
of analyses of the processes it can be stated that the characteri- 
stic frequencies of the processes associated with the typical 
use of the engine are contained within the interval (0.01 + 10) 
Hz, with a large safety margin [3]. 

The slowly changeable processes in comparison with those 
connected with the typical use of the engine are of entirely 
different dynamical character. To this category can be counted 
the thermal processes which characterize engine’s thermal sta- 
te described by the set of temperatures of engine parts and sys- 
tems. Such processes have time constants in the order of seve- 
ral or a dozen or so minutes, a few dozen of seconds at least. 
The processes even more slowly changeable than the thermal 
ones are those describing atmospheric conditions , all the more 
tribological processes associated with engine’s wearing [3]. 

In the combustion engine it is characteristic — due to its 
working principle — that the fast changeable processes, espe- 
cially those connected with working cycles, are pseudo-perio- 
dical. The pseudo-periodicity, but not periodicity, of the pro- 
cesses results : firstly from the fact that periodic function is 
determined for unlimited set of independent variable values, 


1) The conventional notion of ,,characteristic frequencies” is assumed equivalent to the range of frequency values for which amplitudes of 
frequency representation of a process have significant values arbitrarily assumed by a researcher (and sometimes customary assumed by 


specialists to be valid) [3]. 
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and secondly from that the investigated processes are real. Due 
to the pseudo-periodicity of the processes associated with par- 
ticular engine working cycles, high-frequency disturbances of 
signals can be modelled as the signal amplitude disturbances 
and also as independent variable ones. In this case the inde- 
pendent variable is usually : 


time - t ; crank angle - a. 
Let x(t) be a pseudo-periodical signal such that : 
|x(t)-x(t+j-T)|<e (1) 
where : 


particular phase j = 1, 2, ...,N 
e>0 
T - a quantity estimated as a tested signal pseudo-period. 


Therefore the inaccuracy £ depends on : 


+ high-frequency noise amplitudes 
+ characteristics of the set of successive pseudo-periods, {T;}, 
where :i= 1, 2,..., N-1. 


The synchronous averaging of the signal x(t) consists 
in determining the function [1, 12, 13]: 


N 
x(t)=—¥ix(t+}-T) 2) 
N j=0 

On the assumption that values of elements of the set {Tj} 
are constant, the synchronous averaging of the signal makes it 
possible — in compliance with the laws of large numbers (in- 
equality of Chebyshev, laws of Markov and Chintsin [2, 9, 17]) 
— to significantly limit high-frequency noise content in useful 
signal, for a certain class of processes describing high-frequency 

noise, e.g. normal ones : 


lim|[x(t)-x(t)]=0 3) 


The most important limitations of effectiveness of the syn- 
chronous averaging of signals are the following problems : 
determination of a tested signal pseudo-period as well as of 
features of the set of particular pseudo-periods. 


In the case of signal digital processing the accuracy 
of the pseudo-period estimation depends on [1, 3, 4, 12, 13]: 


> resolution of quantization of the signal, that determines 
accuracy of assessing independent variable values 

> sampling frequency which determines Nyquist frequency 
of an observed signal 

> signal observation time deciding on frequency resolution 
of an observed signal. 


Influence of sampling frequency on accuracy of pseudo- 
-period estimation is solely theoretical because it is only re- 
quired the sampling frequency to be at least two times greater 
than that corresponding to expected value of pseudo-period. In 
order to obtain a small inaccuracy value of pseudo-period esti- 
mator it is necessary to use a large value of quantization reso- 
lution as well as long observation times [1, 3, 4, 12, 13]. 

To satisfy the conditions is sometimes really difficult, e.g. 
when the conditions not directly accounted for in testing pro- 
gram influence qualification of processes as stationary ones — 
— in this case an increase of observation time can make it im- 
possible to qualify a process as stationary one. 

The other effectiveness limitation of the synchronous ave- 
raging of signals depends on how large deviations of particular 
pseudo-periods from a value estimated as representative for 
a tested signal, are. In the case of too large deviations the syn- 
chronous averaging of signal can lead to a significant limita- 


tion of amount of information about useful signal, hence — in 
fact — to generation of additional noise instead of limitation of 
noise content in a tested signal (it is obvious that any proces- 
sing of signal inevitably contributes to noise generation [1, 3, 
4, 12, 13]). 

Effective improvement of quality of the synchronous ave- 
raging of pseudo-periodical signals occurring in piston com- 
bustion engine can be obtained by signal sampling in the do- 
main of crankshaft rotation angle, but not in the domain of 
time [1, 3, 4, 12, 13]. This way influence of rotational speed 
fluctuations on effectiveness of the synchronous averaging of 
signals may be avoided. 


INVESTIGATIONS OF THE PROCESSES 
OF INDICATED PRESSURE 
AND INJECTION PRESSURE 


The presented basic problems associated with research on 
fast changeable processes occurring in diesel engines are here 
exemplified by the processes of indicated pressure and fuel 
injection pressure in cylinder. The investigations have been 
carried out at the laboratory stand equipped with 6AL 20/24 
Sulzer ship diesel engine of 37.7 dm” piston displacement, 
420 kW rated output and 750 min ‘rotational speed. 

The tests were carried out in statical working conditions, 
i.e. those independent on time [3] in the frequency range corre- 
sponding to real engine working conditions. In Fig.1 the test- 
ing point of the engine in question is presented. 
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Fig. 1. Testing point of the tested engine operating in compliance with 
rotational speed characteristics 


For testing the signals quantized by means ofa 12-bit trans- 
ducer and recorded at the sampling interval At = 50 Us , the 
sample sets of 32 000 elements each were used. 


The analyses were performed 
in the domains of [1, 3, 12, 13, 17]: 


> time 
> process values 
> frequency. 


The recorded signals are presented in Fig.2 and 3. 
10 
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Fig. 2. Indicated pressure pg in function of time 
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Fig. 3. Fuel injection pressure pw in function of time 

In Fig. 4 the correlation relationship of the indicated pres- 
sure and injection pressure is presented. The correlation coef- 
ficients are of relatively large values : beginning from 0.2441 
of Kendall’s t and 0.2471 of Kendall’s y [2,9], 0.2855 of Spe- 
arman’s range coefficient [2,9] up to 0.6079 of Pearson’s line- 
ar coefficient [2, 9, 14], for this reason there is no ground to 
reject hypotheses on lack of correlation between analyzed sig- 
nals (Fig.5). 
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Fig. 4. Correlation relationship of injection pressure pw 
and indicated pressure Pg 
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Fig. 5. The correlation coefficients of Spearman : R , Kendall : Y and t, 
and Pearson : p - for sets of indicated pressure and injection pressure, 
as well as the probability of non-rejection of hypotheses 
on lack of correlation of the sets , P 


Spearman 


The tested signals were subjected to stationarity analysis 
[1, 12, 13, 17]. To this end the mean values and standard devia- 
tions were determined in particular phases corresponding to 
pseudo-periods (Fig.6 and 7). On the basis of the performed 
analysis the assumption on stationarity of the tested signals 
was adopted. 

The signals of indicated and injection pressures were sub- 
jected to time and frequency analyses. For the investigations 
an algorithm based on Fourier’s fast transformation as well as 
Hamming’s five-sample-wide window was applied [1, 12, 13]. 
The signal power spectral density of indicated pressure and 
injection pressure, with eliminated linear trends, are presented 
in Fig.8. 
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Fig. 6. Mean value AV and standard deviation D 
of indicated pressure pg in course particular phases 
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Fig. 7. Mean value AV and standard deviation D 
of injection pressure pw in course particular phases 
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Fig. 8. Signal power spectral density G 
of indicated pressure pg and injection pressure Pw 


The first maximum of power spectral density of both signals 
appears at the frequency f = 5.4931640625 Hz, which corre- 
sponds to the pseudo-period T = 1/f= 0.182044 s. The pseudo- 
-period of that value corresponds to the rotational speed : 


n =f -2-60 [min] (4) 
where : frequency f [Hz] 


On substititions : n = 659.2 min | the mean rotational spe- 
ed recorded during the test series at the frequency of 1 Hz, was 
equal to 659.3 min — [6]. On the basis of the results, inaccura- 
cy of pseudo-period estimation of the tested signals can be as- 
sessed : it is smaller than 0.023%. 


Also, the auto-correlation function of indicated 
and injection pressures was determined (Fig. 9). 


Also, the results of the correlation analysis of signals con- 
firmed the determined pseudo-period value : the maximum of 
correlation function occurred at the time-lag At = 0.182044 s. 
The mutual testing of the signals of indicated and injection 
pressures was also performed (Fig.10 and 11). 


*) The injection pressure Pw was measured in voltage units and due to technical reasons it was not converted into pressure units 
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Fig. 11. Mutual correlation function Rm of indicated Pg 
and injection Pw pressures 


This mutual testing of the signals confirmed the estimated 
pseudo-period value representative for the tested signals. 

By making use of the determined pseudo-period value of 
the tested signals their synchronous averaging was performed. 
Before the synchronous averaging the courses were subjected 
to the lowpass filtering of 1 kHz limit frequency to decrease 
high-frequency noise content in the tested signals. For realiza- 
tion of the operation a filter with applied Fourier’s fast trans- 
formation algorithm, was used [1, 4, 5, 12, 13]. 

In Fig.12 a fragment of the initial course and the averaged 
one of indicated pressure, and in Fig. 13 — of injection pressure, 
was presented. 

In both the cases the mean deviations of initial courses are 
rather small. The courses of the deviations in particular phases 
are presented in Fig.14 and 15, and the sets of all deviations — 
— in Fig.16 and 17. 

Despite the relatively small deviations of signal values from 
the averaged one, the interference phenomenon is here charac- 
teristic, especially distinct for injection pressure the course of 
which shows content of high-frequency components of a gre- 
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Fig. 12. A fragment of the initial course 
and the averaged one of indicated pressure Pg 
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Fig.13. A fragment of the initial course 
and the averaged one of injection pressure Pw 
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Fig. 14. Fragments of the courses of deviations from the averaged value 
of indicated pressure Apg in particular phases 
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Fig. 5. Fragments of the courses of deviations from the averaged value 
of injection pressure Apw in particular phases 


ater amplitude than in the case of indicated pressure. The inter- 
ferences occuring in the sets of deviations are due to : on one 
hand — real non-periodicity of the tested courses, on the other 
hand — very small difference of estimated pseudo-period and 
successive real pseudo-periods. 
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Fig. 16. Set of deviations from the averaged value 
of indicated pressure Apg in particular phases 
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Fig. 17. Set of deviations from the averaged value 
of injection pressure Apw in particular phases 
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In Fig.18 and 19 are presented the results of analysis of 
standard deviation of the tested signals and deviations from 
their averaged value. Small values of standard deviations and 
deviations of these deviations, relative to the standard devia- 
tions of the initial signals, were observed : smaller — of about 
1.5% - for indicated pressure, and greater — of about 5% - for 


0.02176 


D(pg) , D(Apg) [MPa] 
D(pw) , D(Apw) [V] 


0.000328 0.000197 


D(Pg) D(Pw) D(Apg) DAP w) 
Fig. 18. Standard deviation D of the sets : of indicated pressure Pg 
and of injection pressure pw , as well as of deviations 
from averaged values of these pressures Apg , Apw 


D(Apg)/D(g) , D(Apw)/D(pw) [-] 


Pg Pw 
Fig. 19. Ratios of standard deviations D of the sets: of deviations 
from averaged values of indicated and injection pressures, (Apg , Apw ) 
as well as of indicated pressure pg and injection pressure Pw 
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injection pressure — in this case a greater content of high- 
-frequency noise not eliminated from the initial signal, can be 
clearly observed (Fig. 20). 
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Fig. 20. Power spectral density G of standardized deviations 
from averaged values of indicated and injection pressures (Apgs , Apws ) 


Also, correlation relationship of deviations of the standar- 
dized indicated pressure (of the mean value equal to zero and 
standard deviation equal to | [1, 2, 9, 12, 13]) and deviations 
of the standardized injection pressure, was determined. The 
correlation coefficients of Spearman, Kendall and Pearson of 
the analyzed sets , as well as the probability of non-rejection of 
hypotheses on lack of correlation of the sets are presented in 
Fig.21 [2, 9, 14]. 


Pearson 


Gamma Tau 


Spearman 


Fig. 21. Correlation coefficients of : Spearman, R, Kendall, y and T, as well 
as of Pearson, p, of the sets of standardized deviations from averaged 
values of indicated pressure Apgg and injection pressure Apws , 
and the probability of non-rejection of hypotheses on lack 
of correlation of the sets, P 


On the basis of the performed analysis it was stated that 
there is no ground to reject hypothesis on lack of correlation 
between the sets of standardized deviations from averaged va- 
lues of indicated pressure and injection one. 

The sets of deviations of the pressures in question as well 
as standardized deviations were analyzed in the value domain. 
In Fig.22 the histograms of standardized sets of deviations of 
indicated pressure and injection one from their averaged valu- 
es, are presented, and in Fig.23 — the probability density of the 
tested sets. 

Also, conformity of the tested sets to normal distribution 
was assessed. To this end the following hypotheses were ap- 
plied: of Kolmogorov and Smirnov [2, 9, 10, 16], of Lilliefors 
[2, 9, 11], and of Shapiro and Wilk [2, 9, 15]. In Fig.24 for 
instance the Kolmogorov—Smirnov statistics of the sets of the 
standardized deviations of indicated pressure and injection pres- 
sure from their averaged values, are presented; the probability 
of non-rejection of the above mentioned hypotheses equals zero. 

As a result of the tests of the analyzed sets it was stated that 
there has been no ground to accept the used hypotheses on 
conformity to normal distribution of the sets of standardized 
deviations of indicated and injection pressures from their ave- 
raged values. Probably, a big interference content in the analy- 
zed sets to a large extent decided on the observed, different 
from normal, character of the analyzed deviations. In the case 
of strictly synchronous averaging the tested signals, the nor- 
mally distributed deviations from averaged value should be 
rather expected. 
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Fig. 22. The histograms of sets of the standardized deviations of indicated 


pressure , Apgs , and of injection pressure, Apws , from their averaged 
values (continuous lines stand for respective normal distributions) 
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Fig. 23. The probability density g of sets of the standardized deviations 
of indicated pressure, Apgg , and of injection pressure, Apws , 
from their averaged values 
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Rys. 24. The Kolmogorov—Smirnov statistics d of the sets 
of the standardized deviations of indicated pressure Apg 
and injection pressure Apy from their averaged values 
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SUMMARY 


Investigations of fast changeable processes make great dif- 
ficulties which can result in over-interpretation of test re- 
sults, namely in formulating judgements on investigated 
phenomena, on the basis of the investigation of features of 
measurement systems, as well as on the knowledge gained 
during analyzing the phenomena not directly accounted for 
in the test program. Therefore a formal critical analysis of 
measurement results, based on real knowledge on investi- 
gated phenomena as well as formal tools for processing the 
measurement results, is necessary. 


The example analyses presented in this paper have confir- 
med that it is reasonable to comprehensively test recorded 
signals as well as noise features. Such tests can a.o. to justi- 
fy the thesis on a physical nature of noise eliminated from 
the initial signal. As a result of the analysis it can be also 
stated that the transformations performed on a recorded sig- 
nal may be a factor deciding on the nature of noise. Such 
case was also revealed in this work: some features of the 
tested noise originated from the process of selecting it out 
of the initial signal, namely the interferences due to syn- 
chronous averaging. 


The comprehensive testing of recorded signals is necessa- 
ry to avoid subject-matter errors difficult to be predicted in 
further analyses. Examples of such errors resulting from 
transformations of the signals containing noise of an un- 
known nature are numerous in the field of combustion en- 
gines; for instance, the determination of heat rate generated 
in engine cylinder, performed on the basis of indicated pres- 
sure signal. In this case noise due to numerical differentia- 
tion of indicated pressure is usually subjected to such pro- 
cedures [4, 5]. 


NOMENCLATURE 
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— mean value 

— Kolmogorov—Smirnov statistics 

— standard deviation 

— frequency 

— probability density 

— spectral density of power 

— mutual spectral density of power 

— occurence frequency of a process value on histogram 
— effective torque 

— rotational speed 

— indicated pressure 

— injection pressure 

— probability of non-rejection of a hypothesis 
— Spearman’s range correlation coefficient 
— auto-correlation function 

— mutual correlation function 

— time 

— signal pseudo-period 

— process 

— synchonously averaged process 

— time-lag 


— indicated pressure deviation 
— standardized deviation of indicated pressure from its 
averaged value 


Apy — injection pressure deviation 


Apws — Standardized deviation of injection pressure from its 
averaged value 

Yy — Kendall’s gamma correlation coefficient 

p — Pearson’s linear correlation coefficient 

T — Kendall’s tau correlation coefficient 
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C onference —* 
COMPOWER'04 


On 2-3 December 2004 Gdansk University of Tech- 
nology and Alstom Power Ltd organized the 3rd Interna- 
tional Scientific Symposium on : 


Technical, Economic and Environmental 
Aspects of Combined Cycle Power Plants 


It was first of all aimed at drawing attention of engine- 
ers and economists onto application of turbines to propul- 
sion and heating systems, especially those based on so 
called gas-steam cycles. They are characterized by tech- 
nical, economical and ecological advantages such as : 


> the improvement of the power plant performance 
> the increase of the energy conversion efficiency 
> the decrease of the coal consumption 

> the decrease of the environmental pollution. 


29 papers, split into 5 topical chapters, 
were prepared for presentation during the Symposium : 


> Modelling and simulation of turboset characteristics 
(6 papers) 

> Neural network application to turboset simulation 
(6 papers) 

> Turbine plant diagnostics (4 papers) 

> Analysis and improvement of turbine plant characte- 
ristics (8 papers) 

> Turbine plant operation, reliability and control 
(5 papers). 


The papers were mainly prepared by Polish scientific 
workers from : Institute of Fluid-Flow Machinery of Po- 
lish Academy of Sciences, Gdansk University of Techno- 
logy, Silesian University of Technology, Warsaw Univer- 
sity of Technology, Technical University of Zielona Gora, 
Gdynia Maritime University, Naval University of Gdy- 
nia. 3 papers were elaborated by foreign authors : from 
Mexico, Russia and Sweden. Employees of Alstom Po- 
wer Ltd much contributed to the Symposium, and presen- 
ce of representatives of two Polish firms : Polteknik and 
Diagnostyka Maszyn could be also clearly observed. 


An interesting part of the Symposium 
was Round Table Discussion on : 


"Turbine Industry and Energetics Dilemmas 
(efficiency, reliability, diagnostics, control)" 


A very attractive accent of the scientific meeting was 
the trip to Władysławowo, a Baltic port, where the Sym- 
posium's participants had the occasion to visit a gas heat- 
-generating plant. The modern object built in 2003 takes 
natural gas, through the longest-in-Europe underwater pi- 
peline, from a gas production rig 82 km distant offshore. 


MARINE ENGINEERING 


Prediction of corrosion fatigue crack 
propagation life for welded joints 
under cathodic potentials 


Marek Jakubowski 
Gdansk University of Technology 


ABSTRACT 


Enhancement of corrosion fatigue crack growth rates by cathodic protection is observed 
below the optimum applied potential of the protection. An empirical formula for the 
effect of the protective potential below -0.8 V (Ag/AgCl) on the crack growth rates for 
some classes of shipbuilding steels tested in salt water has been derived for medium and 
high value ranges of the stress intensity factor. For the lower value range the formula 
reflects a relatively steep decrease of the crack growth rates (against the same values in 
air) along with decreasing the stress intensity factor range. A simple formula for the 


corrosion fatigue crack propagation life under cathodic potential has been derived for fillet welded 
joints in bending by integrating the corrosion fatigue kinetic characteristics. The new formula for the 
stress intensity factor range has been used for the bent joint. The predicted "S-Np" curves have been 
compared with experimental data, taken from literature, for two different values of both applied potential 
and the plate thickness. The predicted curves correspond approximately to lower bounds of the test 
results. The presented procedure can be applied to joints of higher strength steels (of the yield strength oy 
= 315 and 355 MPa) fatigued at any applied cathodic potential below -0.8 V under sea loading 
of (0.05+0.2 Hz) frequency at (0+0.2) stress ratios. 


Keywords : corrosion fatigue, shipbuilding steel joints, cathodic protection 


INTRODUCTION 


One of the main aims of fatigue investigations is to create 
fatigue life calculation procedures for structural members. In 
the paper [1] this author presented a simple procedure for jo- 
ints fatigued in air and verified it by using fatigue test results 
for notched steel specimens. The specimen geometry in the 
notch region simulated that of fillet welded joint. In the au- 
thor’s paper [2] the "S-N," curves calculated by using the pro- 
cedure were compared with the fatigue test results for real fil- 
let welded joints. The procedure gave rather conservative pre- 
diction of the total fatigue life, except for the largest values of 
the stress concentration factor K; , and the longest life region 
(above 1+3x10 cycles), hence it can be recommended for ap- 
proximate fatigue life predictions. However, real fatigue life 
of ship and offshore structures is usually affected by marine 
environment. A prediction procedure of corrosion fatigue crack 
initiation life and total life for welded joints under free corro- 
sion conditions was presented elsewhere [3]. 

There is a common opinion that cathodic protection incre- 
ases corrosion fatigue strength and restores true fatigue limit, 
although the real magnitude of this effect is not certain. Some 
tests of notched steel specimens [4 + 6] have shown that the 
strength in salt water under cathodic protection is higher than 


in air, however, other tests [7 , 8] have shown that the strength 
is lower than in air. Therefore the "S-N," curve in air cannot 
be always considered as conservative one for cathodically pro- 
tected structures in salt water. Hartt and Hooper [5 , 6] have 
asserted that the fatigue limit for the protected specimens is 
controlled by hydrogen embrittlement and the cathodic-depo- 
sit-induced crack closure effect. Hydrogen embrittlement is well 
known as a crack-growth-rate accelerating process. Crack clo- 
sure phenomenon leads to reduction of the crack growth rate 
by reduction of an effective value of the stress intensity factor 
that is considered to be the true crack-driving parameter. Thus 
they are competitive processes. However, the deposits forma- 
tion process within the crack and their interaction with the crack 
walls is complex and practically unpredictable. The kind of 
the mentioned crack closure affects especially propagation of 
short cracks that is usually considered as a dominating part of 
a conventional crack initiation period. It makes the prediction 
of the crack initiation life for cathodically protected joints prac- 
tically impossible at present. Therefore this paper deals only 
with the prediction of the corrosion fatigue crack propagation 
life Np for welded joints under cathodic protection. The pre- 
dicted "S-N," curves were compared with experimental data 
evaluated by Vosikovsky at al. [8] for the welded fillet joints 
schematically shown in Fig.1. 
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Fig. 1. Geometry of the analyzed welded joints 


CORROSION FATIGUE CRACK 
GROWTH RATE CHARACTERISTICS 


Fig.2 shows the effect of the applied cathodic potential ® 
on the ratio Beat (at AK = 20 + 40 MPaVm) of the crack growth 
rate (da/dN)@ at the applied potential ®, and the crack growth 
rate (da/dN),i, in air for BS4360-50D steel (oy=360 MPa) [9], 
as well as for St41U5 steel (oy = 316 MPa) and St41E-TF32 
steel (Oy = 318 MPa)[10, 11]. 
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Applied potential, ®, V(Ag/AgCl) 
Fig. 2. The ratio Beat = (da/dN)@/(da/dN)air versus the applied potential ® 


for BS4360-S0D steel, acc.to [9], and for St41U5 steel and St41E-TF 32 
steel, acc. to [10,11] 


For ® < - 0.8 V (Ag/AgCl) the crack growth rate (da/dN)@ 
can be approximated by the following formula [10] : 


da da da 
—— |= —— | =A|——_| - (Œre -®) 0) 
(an, i a(x), (x) n 


air 
where : 
Efree - free corrosion potential 
A -coefficient of the applied cathodic potential effect. 


A equals 12.2 v`! for the upper envelope of the results, 
and 7 V for their lower envelope; and the mean value — about 
9.6 V . Approximate values of A, evaluated by this author for 
a higher strength steel (oy=370 MPa) at AK=25+35 MPaVm 
on the basis of the test results [12] are 3 + 14 V`, i.e. they are 
of the same order as the above mentioned values. 

Some tests results [9 + 14] showed that the strongest en- 
hancement of the crack growth rate by hydrogen embrittlement 
for steels under cathodic protection occurs approximately for 
AK above 20 MPavVm. Another set of data [7] leads to the same 
conclusion. The above is true for the loading frequencies f = 
= 0.1 + 0.2 Hz and the stress ratios R = 0 + 0.2. Therefore it is 
assumed that the crack growth rate at any applied potential 
® <- 0.8 V and for AK > 20 MPavm is given by the equation 
(1) with A= 9.6 Vv. For AK < 20 MPaVm hydrogen embrit- 
tlement is not so dangerous as for AK values greater than 
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20 MPavm, and in this case the crack growth rate quickly drops 
to the values observed in air, or below them. It is ad hoc assu- 
med that (da/dN),i=(da/dN)@ for AK = AKo . Value of AKo 
corresponds to the end of conventional crack initiation stage 
and to the beginning of conventional propagation stage [1 + 3]. 
Fig.3 schematically shows the assumed corrosion fatigue crack 
propagation characteristics. 


da/dN 


Fig. 3. Schematic diagram of the assumed corrosion 
fatigue crack growth characteristics under cathodic potentials 


STRESS INTENSITY FACTOR 


The stress intensity factor range AK 
was evaluated by using the formula : 


AK = AS./na- Y(a/t,) (2) 
In the present calculations, and in [1 , 2] as well, the follo- 
wing form of the correction function Y was assumed : 
-Ô 
Y =y; (a/tw) (3) 


For smooth beam in bending : Y = 1.12 for a/ty = 0 and 
a/ty = 0.305. For notched beam the following values were as- 
sumed : Y(0.000305) = 1.12K; and Y(0.305) = 1.12, that led 
to the following equations : 


log(K,) 
ae $ 
y= 1.12 (0.305) (5) 


The stress concentration factor for the joint shown in Fig.1 
was calculated by the following formula [15] : 


K; = 1 + 0.512 0 °°”? (t,/1)°40? (6) 


The values of Y calculated by using the above given for- 
mulae are comparable [2] to those evaluated by Niu and Glin- 
ka [15] by means of the finite element method. 


THE ASSUMED CHARACTERISTIC 
CRACK LENGTHS 


Values of aj are usually — by a convention — assumed equal 
to 0.5 + Imm. In this paper, as well as in [1 , 2], a; values — by 
another convention — correspond to a fixed value of AK =AKg , 
where AKọ corresponds — by an ad hoc assumption — to the 
end of the near threshold propagation rate (da/dN=5x107'm/c). 
Values of aj were calculated by means of the formula : 


AK, 


1/0.5-8 
i -( TPY PAS A 


ai < 0.0014 m (7b) 

The final crack length is assumed ar = 0.4ty, , i.e. the same 

as in [2] for similar welded joints fatigued in air. The value of 

ar was evaluated from the intersection of the adjacent seg- 

ments of the corrosion fatigue crack growth rate characteri- 
stics (Fig.3), i.e. by the following equation : 
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CORROSION FATIGUE CRACK 
PROPAGATION LIFE 


Integration of the assumed crack growth rate characteristics 
(Fig.3) leads to the following formula : 
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where : 


a = 1- m, (0.5 - 8) 
o =1- ee 5) 
2 =5°10 
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Fig. 4. Comparison of the corrosion fatigue crack propagation curves, 
"S-Np" , acc. to equation (9) with the tests results published in [8] 
(here shown as the points) 


In Fig.4 corrosion fatigue lives of the welded joints shown 
in Fig.1, predicted by the above given formula, are compared 


to the lives of the same joints, determined empirically by Vosi- 
kovsky at al. [8] for two values of the applied cathodic poten- 
tial. In order to evaluate the crack propagation life N, obtained 
by Vosikovsky the crack initiation life N; was subtracted from 
the total life, both read off from the appropriate figures in [8]. 
Values of Np , especially those for ty = 26 mm, demonstrate 
rather large scatter ; it can be partly due to the reading-off er- 
rors. The predicted "S-N," curves give a conservative asses- 
sment as they approximately correspond to the lower bounds 
of the test results, although for thicker joints such prediction is 
only slightly conservative. 


SUMMARY 


* A procedure for conservative prediction of the corrosion 
fatigue crack propagation life N, of cathodically protected 
welded steel joints was presented. The procedure has limi- 
ted applicability — it can be applied to joints manufactured 
of higher strength steels (of class E315 and E355, denoted 
by IACS as AH32, DH32, EH32, AH36, DH36, and EH36) 
fatigued at any applied potential lower than - 0.8 V (Ag/ 
/AgCl,) under sea wave loading of frequencies within the 
range of 0.05 + 0.2 Hz and stress ratios R = 0 + 0.2. 


* The predictions for joints of steels of lower strength classes 
at any higher loading frequency and R values are expected 
to be more conservative, whereas for steels of higher clas- 
ses (E390 or higher) at a lower frequency and R values it 
should be expected that the procedure can overestimate the 
crack propagation life of the joints. 


* Residual welding stresses were not considered. 


* The final expression for N, concerns fillet welded joints in 
bending, however, the procedure can be applied to every 
welded joint geometry for which it is possible to calculate 
the stress intensity and the stress concentration factors. 


NOMENCLATURE 

a — crack length 

ar — final crack length, [m/cycle = m/c] 

a; — initial crack length, [m/c] 

A — coefficient of the applied cathodic potential effect 

c — cycle 

C — material constant in the Paris equation, [m - (MPa Vm)"] 
Ca — material constant in the Paris equation for region I under 


cathodic protection, [m : (MPa Vm)" `] 


daj fatigue crack growth rate in air, [m/c] 
dN Jair 
OF ja corrosion fatigue crack growth rate under cathodic 
dN © potential ®, [m/c] 
Efee — free corrosion potential, [V](Ag/AgCl,) 
AK  — range of stress intensity factor , [MPa Vm] 
AK; — final stress intensity factor, [MPa Vm] 
Ki — stress concentration factor, [-] 
m — exponent in the Paris equation in air, [-] 
ma — exponent in the Paris equation for region I under cathodic 
protection, [-] 
N — fatigue life, [cycles] 
Np — crack propagation life, [cycles] 
P — loading force, [kN] 
r — weld toe radius, [m] 
R — stress ratio (= Syin/Smax), [-] 
S — nominal stress, [MPa] 
AS — nominal stress range, [MPa] 
tp, tw — thickness of joined plates, [m] 
Y — correction factor in the formula for stress intensity factor, [-] 
Q, &, — functions appearing in eq.(9), [-] 
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Beat — coefficient of cathodic potential effect on crack growth 


rate, [-] 
Y — constant in eq. (3) [-] 
ô — exponent in eq.(3) [-] 
Oy — yield point stress, [MPa] 
p — cathodic potential, [V] (Ag/AgCl,) 
8 — weld toe angle, [radian] 
IACS — International Association of Classification Societies 
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POMERANIAN SCIENTIFIC CONFERENCE 


During the Conference were presented and discussed 
49 papers prepared by specialists from 11 scientific rese- 
arch centres including Hungarian Academy of Sciences, 
Budapest. 


To the Conference program 
most contributed scientific workers of : 


¢ Gdansk University of Technology 
(22 papers of single authors + 9 of co-authors) 


è Gdynia Maritime University (4 + 4) 

è Koszalin Technical University (5 + 1) 

¢ Ship Design & Research Centre, Gdansk (2 + 2) 
Technical University of Szczecin (2 + 0) 

Institute of Fluid-Flow Machinery, Gdansk (0 + 3). 
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Miscellanea 


Combustion in Piston Engines 


Technology « Evolution « Diagnoses « Control 


A very inspiring book written by Prof. Antoni K. Op- 
penheim of University of California, has been recently pu- 
blished by Springer-Verlag Berlin Heidelberg New York 
(marked ISBN 3 - 540 - 20104 - 1). 


As the title of the book suggests it concerns the techno- 
logy of combustion in internal combustion engines - the pri- 
me automotive powerplant today. The problem is important 
because of the escalating emission of pollutants by combus- 
tion and the dwindling natural resources of energy. As 
a consequence, the automotive industry found itself confron- 
ted with a significant amount of societal pressures and go- 
vernmental regulations, which brought about a state of flux 
and confusion. 


Confronted with the mounting dissatisfaction with cur- 
rent technology of internal combustion engines, the auto- 
motive industry is involved nowadays in a frantic search for 
alternative fuel and alternative energy conversion systems. 
The question how much effort should be spent on searching 
for an alternative system, while the attractive prospects of 
equivalent progress in the technology of combustion in pis- 
ton engines are disregarded, is open. 


The purpose of the book — its author says in Preface — is 
to provide the background for implementation of such a pro- 
gress. This is accomplished by providing an assessment of 
the technology of combustion in piston engines, followed 
by furnishing an engineering method of approach for evalu- 
ation of the effectiveness with which fuel is utilized in the 
engine cylinder, and demonstrating how it can be improved. 


These subjects are exposed, respectively, in two parts : 
Part 1 made of : 


Ch. 1.: Overview, describing what the technology of com- 
bustion in piston engines is all about 

Ch. 2.: Perspective, providing an account of how did it 
got to its present state 

Ch. 3.: Prospective, pointing out how it technology can 


be advanced and what gains can be thereby derived 


Part 2 is based 
on the recognition that combustion 
is at the heart of a piston engine (hence chapters 
of this part are entitled in medical terms) 


Ch.4: Diagnosis introduces the principles of a "engine 
cardiology" profession 

Ch. 5.: Procedure prescribes the manner in which it is ap- 
plied to a given engine 

Ch. 6.: Prognosis provides an insight into its prospects. 


In the author's opinion the technology of combustion in 
piston engines is today at the threshold of a quantum leap, 


comparable to that which took place in electronics by its 
transition from vacuum tubes to transistors, and that his 
book provides a recipe for advancement of the technology 
in question. 


Engines of the future, based on this advancements, wo- 
uld feature distributed combustion executed by means of 
Pulse Jet Injection and Ignition ( PJI&I) actuators modula- 
ted by Micro-Electronically Controlled Combustion Systems 
(MECC) in response to signals provided by pressure sen- 
sors to assess the effectiveness with which fuel is utilized in 
an engine cylinder, and to monitor the operation of MECC 
system in optimizing the performance of a PJI&I system. 


The book suggests that by the use of these systems it 
would be possible to build a car engine for which the forma- 
tion of pollutants can be reduced by orders of magnitude, 
while the consumption of fuel is halved. 


The book, written in a very clear and concise way, con- 
sists of 160 pages and includes many figures and diagrams 
and a dozen or so reference sources. 


Combustion 
in Piston Engines 


© Springer 
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OPERATION & ECONOMY 


Application of statistical methods 
and artificial neural networks 
for approximating ship’s roll in beam waves 


Tomasz Cepowski 
Maritime University of Szczecin 


ABSTRACT 


In the paper, were presented approximations of numerical calculations of ship ’ roll on the 
basis of main service parameters of the ship. This way were obtained several relationships 
which make it possible to approximate ship roll in regular and irregular waves by using 


<7 = 


the parameters available in the phase of voyage routing. The relationships were elabora- 
ted by means of artificial neural networks as well as linear and non-linear regression 
>, methods. A comparative analysis of the methods regarding approximation accuracy 


against standard data was also performed. 


Keywords : ship's roll, regular wave, irregular wave, approximation, neural networks, regression methods 


INTRODUCTION 


Within the contemporary shipping problems there is a pro- 
blem of searching for ship voyage route satisfying several crite- 
ria among which the following usually are the most important : 


> voyage duration time (at assumed parameters of ship pro- 
pulsion system) 
> operational cost which mainly depends on fuel consumption. 


A limitation which seriously influence route choice is the 
condition of safe shipping, which consists of many factors 
among which ship’s behaviour in rough seas should be distin- 
guished. 

So many parameters influence ship’s behaviour in waves 
that accounting for all of them makes the voyage routing pro- 
cess very complicated. Hence, out of all the parameters, ship 
master is forced to take into account only most important ones, 
among which wave parameters, ship’s motion parameters or 
also selected hydromechanical parameters of the ship, are num- 
bered. The small amount of available information does not make 
it possible to use exact methods of determination of ship mo- 
tion in waves. 

In the subject-matter literature, methods making it possible 
to solve the problem in a satisfactory way, are still lacking. 
Design recommendations dealing with ship seakeeping quali- 
ties, given in classification rules, are of a very limited charac- 
ter. The calculation procedures there presented are rather inac- 
curate and they allow to determine only ,,designed” ship mo- 
tion amplitudes and accelerations which are connected only to 
a certain degree with real ship behaviour in waves. The appro- 
ximations of ship seakeeping qualities published in the scienti- 
fic literature are too general, rather inaccurate, and usually ap- 
plicable only to a given hull form [1,5]. 
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In this paper an attempt to solve the problem has been un- 
dertaken, aimed at elaboration of a simplified but exact model 
of predicting ship’s roll on the basis of main ship service para- 
meters. 


METHOD 


The research in question was limited to approximation of 
ship’s roll in regular and irregular waves coming from the di- 
rection perpendicular to ship’s plane of symmetry. It was assu- 
med that approximations of the oscillations have to be elabora- 
ted on the basis of the parameters taken into account in ship 
voyage routing and simultaneously having significant influen- 
ce on the motions. Among such parameters the following are 
usually numbered : 


+ wave parameters in the form of its height and period (in the 
case of aproximation of ship roll in irregular waves) 

+ ship motion parameters 

+ ship service parameters associated with, a.o., ship mass dis- 
tribution and its hull hydrostatic characteristics. 


In the research in question it was adopted an approach of 
approximating the standard values of the roll angles ¢ (deter- 
mined by means of exact numerical methods) by using the ap- 
proximating function fon the basis of the n-element set of in- 
put parameters W[X), Xo, ... Xp]: 


b = f(W) (1) 
where : 
W - n-element set of input parameters : X4, X2, ... Xp 
þ - standard values of roll angles calculated by me- 
ans of the exact method 
approximating function searched for. 


= 
1 


STANDARD VALUES OF ROLL ANGLES 


Standard values of roll angles were determined by using 
exact numerical methods based on the two-dimensional flow 
theory, i.e. the SEAWAY software (was elaborated by Shiphy- 
dromechanics Laboratory, Delft University of Technology, The 
Netherlands) which calculates ship motions in regular and irre- 
gular waves. In Fig.1 are presented the ship roll amplitude cha- 
racteristics calculated by means of the software, and showed 
together with the characteristics obtained from model tests [10]. 

5 


RAO roll [-] 


ML 
+ Fn=0.1, Exp. Zhou [10] —— Fn=0.1, Seaway 
© Fn=0.22, Exp. Zhou [10] . . - Fn=0.22, Seaway 


Fig. 1. Roll angle amplitude characteristics of the containership 
„Nedlloyd Dejima” of : L = 270 m, B = 32.20 m, d = 10.85 m, Cg = 0.596, 
ZG = 16.45 m, in beam waves, acc. [10] 


The approximations and numerical calculations were car- 
ried out for the S-175 model containership of the following 
main dimensions : 


L (length between perpendiculars) = 175 m 
B (breadth) = 25.4 m 
d (design draught) = 9.5 m. 


Its hull form is given in Fig.2. 


Fig. 2. Theoretical hull frames of S-175 containership, acc. [4] 


The numerical calculations were performed by applying, 
a. o., the modified method of two-dimensional flow, accoun- 
ting for diffraction of wave loads [4] , and assuming ITTC wave 
spectrum. 

On the basis of literature data, were selected the parame- 
ters which have decisive influence on ship roll, and simultane- 
ously considered as the service ones. To them it belongs, a.o. : 


= the ship volumetric displacement V 
© the initial transverse metacentric height GM 
= the height of the ship gravity centre Zg. 


Taking into account that the volumetric displacement V can 
be expressed in the form of the mean draught d , and that the 
initial transverse metacentric height GM depends on the hei- 
ght of ship gravity centre Zg, one decided to determine the 
ship roll approximations by using : 


* the mean draught d 
* the initial metacentric height GM 


for which the following ranges 
of their values were assumed : 


A ford: from 7 to9m 
A for GM : from 0.3 to 1.5 m 


and 
A for the ship speed V : from 0 to 20 knots. 


Additionally, for approximation 
of ship roll in irregular waves it was assumed : 


# values of the significant wave height Hs : 
from 2 to 4.5 m, taken with 0.5 m step 

Æ values of the characteristic wave period T : 
from 6.5 to 14.5 s, taken with 2 s step. 


For all the above mentioned variants roll angle amplitudes 
of ship in regular and irregular waves were calculated, and this 
way the set of standard data for approximations was obtained. 


APPROXIMATION OF ROLL AMPLITUDES 
OF S-175 CONTAINERSHIP 
IN REGULAR WAVES 


The approximation of ship roll amplitudes in regular wa- 
ves (in irregular waves as well) was carried out by means of 
statistical methods, i.e. linear and non-linear regression, as well 
as artifical neural networks. The latter method belongs to rela- 
tively novel mathematical-numerical methods placed in the field 
of „artificial intelligence”, and it finds more and more applica- 
tions in different domains of science and engineering. The re- 
search results published in [6, 8, 9] show broader and broader 
possiblities of application of artificial neural networks to the 
problems of ship design and operation. 

On the basis of [1,8] and the performed investigations it 
has appeared that the best model of roll amplitude characteri- 
stics against the assumed standard data is the following : 


Yo =] 
(2 (2) 
6) 


where : 


a,b -— coefficients dependent on ship service parameters 
(00) — wave frequency [s]. 


The coefficients a and b were approximated by using the 
statistical methods and artificial neural networks, as shown 
below. 


Approximation of the coefficients a and b 
by using statistical methods 


To approximate the coefficients a and b in the equation (2) 
linear and non-linear regression was used. When using the li- 
near regression the following relationships appeared the best 
out of all investigated models : 
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a=O,d+a,,GM+a,-d* + 
M (3) 
ras bay Vay GM 


2 
b=a,-d+a.- = + Oly: d’ (4) 


where : 


Q0, ..., Qg — coefficients whose values are given in Tab. 1 
The remaining symbols — as explained above. 


Table 1 


fo. 17-09. se 0390035 of 05 


When applying non-linear regression 
the following relationships appeared the best ones : 


œ +0,-d4+ Q, V+ 


5 
tay GM +04," G 


a = Œy +exp 
b=b Bo +b,-B, (6) 
where : 
by = (a, +0; -d + Og: V +0 GM) 
b, = (Œ +,,:d+0,,-V+a,,-GM) 
B, =1 if by <Q, otherwise Ba =0 
B, =1 ifb, >Q, otherwise B, =0 


... Q14 — coefficients whose values 
are given in Tab.2 and 3. 


Qo, 


Table 2 


ee 


Table 3 


so | ea 


sao [oon [oos [nore isa ooz [20 | 


Approximation of the coefficients a and b 
by means of artificial neural networks 


Among all investigated networks used for approximating 
the coefficients a and b in the equation (2) the MLP neural net- 
works of the below given form, showed the highest accuracy : 


l C 
1+¢ (4 VMI: S+P) - A-B) -C70 70] 
a,b = tT (7) 
Q, 
where : 


the symbols : d, GM, V — as explained above. 
S — the vector of normalizing values : [0.5 0.05 0.833] 
P — the vector of displacement values : [- 3.5 0 - 0.250] 


The matrices, vectors and constants 
necessary to determine the coefficient a : 


A — matrix of weight values : 


-0.1304 0.8757 -0.9335 -0.0855 0.9179 -0.0139 0.0338 -0.6867 0.5367 -0.3011 0.8537 -0.1561 -0.9966 
0.4534 0.6205 -0.7704 -0.3399 0.4939 0.8646 -0.7616 0.3134 -0.5045 -0.3016 -0.3004 0.6048 -0.0967 
0.1259 0.0096 0.4784 0.9856 0.6298 0.7879 -0.2093 -0.5050 -0.9623 0.2246 -0.8858 -0.2999 -0.1227 
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B — vector of threshold values : 
[-0.93 -0.764 0.847 -0.201 1.02 0.635 -0.568 -0.181 -0.522 0.926 1.109 -0.516 -0.863] 
C — vector of weight values : 
[1.03 -0.025 -0.528 -0.34 -0.469 -0.64 -0.259 -0.053 0.275 -0.588 -0.810 0.804 -0.082] 
Ql, OL}, Q2 — the coefficients having the values : 
Oo = - 0.323, a = - 0.305, A = 0.496. 


The matrices, vectors and constants 
necessary to determine the coefficient b : 


A — matrix of weight values : 
1.1072 2.6231 5.2094 -4.4498 1.7052 
0.5590 2.2139 -1.6990 1.0816 0.8587 
-1.4717 -2.6282 -2.9165 0.8025 -1.6338 

B — vector of threshold values : 

[-0.515 0.813 -1.374 1.575 -0.769] 

C — vector of weight values : 
[1.562 -1.461 -2.731 -3.637 2.197] 

Q0, OL}, Q2 — the coefficients having the values : 
Qo = - 0.773, Oy = - 0.547, O2 = 3.475. 
Assessment of the approximations 


To assess accuracy of the approximations, 
were used RMS values calculated from the expression : 


(8) 
where : 
RMS -error value 
Yow  — Standard values of ship roll amplitudes 
Yọ | —approximated values of ship roll amplitudes 
n — number of considered variants 
m — number of the values approximated 


on the amplitude characteristics. 
The accuracy analysis was carried out for : 


1) values of the input parameters for which approximations 
were elaborated 

2) the value of the initial metacentric height GM = 2 m, 
being out of the range of the elaborated approximations. 


Values of the RMS errors regarding the elaborated rela- 
tionships, for the above mentioned cases, are given in Tab.4. 
And, in Figs.3 to 5 the elaborated approximations are compa- 


red respective to selected variants. 
Table 4 


MLP neural| Linear |Non-linear 
RMS error i i 
network regression regression 


Tr the input parameters ease) 
Tor GM= 2m (ease 2) 


From Tab.4 and Fig.3 to 5 it results that the approximations 
obtained by using the MLP artificial neural network show the 
greatest accuracy regarding interpolation and extrapolation. 


APPROXIMATION OF ROLL AMPLITUDES 
OF S-175 CONTAINERSHIP 
IN IRREGULAR WAVES 


Approximations of numerically calculated significant va- 
lues of roll angle amplitudes in irregular waves were perfor- 
med on the basis of ship service parameters as well as wave 
parameters. To this end also statistical methods (linear and non- 
-linear regression) and artificial neural networks were applied. 


Explanation to Figs 3+5 : 


exact calculations 
$ by SEAWAY software 


— — — - linear regression 


MLP neural network 
“=----- non-linear regression 


_—-—— — 
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Fig. 3. Values of the ship roll transfer function Yọ, 
at V = 10 kts, d = 9 m, GM = 0.3 m 
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Fig. 4. Values of ship roll transfer function Yọ, 
at V = 20 kts, d=7m, GM = 0.3 m 


0.8 


Fig. 5. Values of ship roll transfer function Yọ, 


at V = 10 kts, d = 8 m, GM = 2.2 m 


Approximation of ship roll amplitudes 
in irregular waves 
by means of statistical methods 


When using linear regression, out of all considered 
relationships, the following one appeared the best : 


Di3 =Q +0 :GM-T.-Hs +a, d- GM: 
-V-T-Hs+a,-d°-GM’ +a, GM? 


where : 
4/3 — significant roll angle amplitude [°] 
Oly , Ly , Q2 , Q3 , Q4 — the coefficients 


whose values are given in Tab.5 
d, GM, V — as explained above 


(9) 


T — characteristic wave period [s] 
Hs — significant wave height [m]. 
Table 5 


a | 


- 0.34231 0.09918 - 0.00021 0.01735 - 1.00743 


And, when applying non-linear regression, the exponential 
relationship of the below given form appeared the best out of 
all considered models : 


a, +a,-d+a,-GM+ 
1/3 =O, +exp i z i (10) 
+a,-Vta,T+a,-Hs 
where : 
Og, Q1, %2, Q3, Q4, Q5, Og — the coefficients 
whose values are given in Tab.6. 
Table 6 


pe >a > |e) a) a) a 


In both relationships the values of the determining coeffi- 
cient prove that the proposed relationships well match up with 
the standard data. From variance analysis it results that both 
equations are essential. Simultaneously the regressions of de- 
pendent variables 61/3 indicate that all components of both equa- 
tions (including their free terms) are essential. 


Approximation of ship roll amplitudes 
in irregular waves 
by means of artificial neural networks 


Out of all considered networks, the following neural 
networks revealed the best approximation features : 


> the MLP network of 5 inputs x 11 hidden neurons 
x | output 

> the RBF (Radial Basic Functions) network of 5 inputs 
x 180 hidden neurons x | output. 


Statistically significant roll angle amplitudes 1/3 approxima- 
ted by means of the MLP neural network can be calculated 
with the use of the following equation : 


l C 
1 + e [60M V.T, Hs]: S+P) : A-B) "UC Ay |7 


aa Gs (11) 
where : 
04/3 , d, GM, V, T, Hs — as explained above 


A — matrix of weight values : 


0.1387 0.3845 -0.3169 -0.1804 -0.3674 -0.6876 0.0606 -0.3752 0.2717 0.0268 9.0807 
0.1743 -3.5447 -2.1067 -1.0170 0.1016 -0.8214 -1.3072 0.1762 3.2725 -0.2252 2.8329 
0.3839 -0.9977 0.2728 0.0992 0.5148 -0.2305 -0.3881 1.5009 -0.2839 0.6204 0.1285 
-0.4993 0.9608 -1.3266 -2.4275 1.0287 -0.2497 -3.1174 -0.1543 3.1220 -0.1918 3.7166 
0.3733 0.9096 0.9033 -0.2308 0.3956 -0.0745 0.4807 -1.1985 0.2697 0.6046 -0.1800 
B — vector of threshold values : 
[0.2133 2.9327 -0.6466 -1.6343 -0.2292 0.7501 -2.5607 -1.9524 2.1205 0.9054 2.9620] 
C — vector of weight values : 
[0.2687 -1.7393 -1.0205 -1.7328 0.3104 0.7946 2.5177 -1.0679 -2.2512 0.7164 2.6939] 
S — vector of values : 
[0.500 0.833 0.050 0.125 0.400] 
P — vector of displacement values : 
[-3.500 -0.250 0.000 -0.813 -0.800] 


Q0, 1, Ol — coefficients having the values : 
Og = - 0.2061, a = - 0.017, a = 0.102. 


To calculate values of the ship roll amplitudes 61/3 appro- 
ximated with the use of the RBF neural network one can apply 
the following generalized equations : 
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4/3 = ([d, GM, V, T, Hs]: A+B)-C (12) 


where : 


A - [5 x 180] matrix of weight values 
B — 180-element vector of constants 
C — 180-element vector of weight values. 


Input values for the equation (12) were normalized to be con- 
tained within the interval : 0...1, and the neurons of hidden layer 
were activated by means of an exponential function. Because of 
too large dimensions of the matrix A and remaining vectors it 
was decided not to present values of elements of the matrices 
and vectors as well as values of normalizing coefficients. 


Values of correlation coefficients of both networks show that 
both proposed models well match up with standard values. 


Assessment of the approximations 


To assess accuracy of the approximations, RMS’ values 
calculated from the following expression, were applied : 


RMS'= (0, = 173 y 


m (13) 


where : 


RMS' — error value 


dy  — Standard values of ship roll amplitudes 
01/3. — approximated values of ship roll amplitudes 
n — number of considered variants. 


The accuracy analysis was carried out for : 


1) values of the input parameters for which approximations 
were elaborated 


2) the value of the initial metacentric height GM = 2 m, and 
that of the significant wave height Hs = 3 + 6 m, being 
behind the range of the elaborated approximations. 


Values of the RMS’ errors associated with the elaborated 
relationships for the above mentioned cases are given in Tab.7. 


Table 7 


ME REH Linear | Non-linear 
RMS’ error neural neural ‘ : 
network | network |"egression | regression 


parameters (case 1) 
for the values behind 
0.69° 5.32° 2.53° 


the range of input 
parameters (case 2) 

From the table it results that in both cases the most accurate 
appear the approximations obtained with the use of the MLP 
artificial neural network. The approximations elaborated by 
means of linear regression are loaded with a small error in the 
case of interpolation, and with rather large one as far as extra- 
polation operations are concerned. 


Addtionally, to more precisely verify the proposed rela- 
tionships the falsification method was applied [7]. To this end, 
such ranges of input parameters were searched for which the 
elaborated approximations appeared the least accurate. Re- 
sults of the investigations are graphically presented in Fig.6+ 
+8 (for interpolation), and in Fig.9 (for extrapolation). The 
diagrams confirm that in both cases the approximation with 
the use of the MLP artificial neural network appears the most 
accurate. 
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Explanation to Figs 6+9 : 
=== — RBF neural network 
linear regression 


exact calculations 
= by SEAWAY software 


—— MLP neural network 


non-linear regression 


10 


0 i z : : : 
0.3 0.5 0.7 0.9 1.1 1.3 
Fig. 6. Significant ship roll amplitudes $173, 
at GM=var,d=7m,V=0,T=14.5s, Hs =4m. 
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Fig. 7. Significant ship roll amplitudes 1/3, 
atd=var, GM=1.5m,V=0,T=14.5s,Hs=4m 
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2 
V [kts] 
0+ r r - 
5 10 15 20 
Fig. 8. Significant ship roll amplitudes 61/3, 
atV=var,d=9m,GM=15m,T=14.5s,Hs=4m 
207 — 


10 15 
Fig. 9. Significant ship roll amplitudes $173, 
atV=var,d=9m,GM=2m,T=95s,Hs=6m 


ACCURACY ASSESSMENT 
OF THE ELABORATED APPROXIMATIONS 
WITH A VIEW OF PREDICTING SHIP ROLL 
IN IRREGULAR WAVES 


The above presented relationships make it possible to ap- 
proximate ship roll in regular and irregular waves. The first 
method allows to approximate roll amplitude characteristics 
on the basis of the ship service parameters (d, GM, V) and next 
to calculate any statistical values of roll amplitudes for any 
wave spectrum. And, the second method makes it possible to 
directly approximate ship roll amplitudes in irregular waves. 

In this part of the investigations the two approaches were 
compared regarding their accuracy in predicting ship roll mo- 
tion in irregular waves against the standard values. The analy- 
sis was performed on assuming the values of input parameters 
to be contained within : 


“* the ranges of the parameters for which the approximations 
were elaborated 

“* behind the ranges of the parameters, in order to test this 
way extrapolation capabilities of the elaborated relationships. 


In the first case the calculations were carried out for : 
d=8m GM=09m Hs=4m. 
And in the second case for : 
d=9.5m GM=2m Hs=6m. 
In both cases the characteristic wave period T was assumed 
equal tol0 s, and the ship speed V within the range from 0 to 
20 knots. 


The analysis was performed for the approximations whose 
accuracy was the greatest and which were described by : 


> for ship roll in regular waves — the relationship (2) having the 
coefficients approximated with the use of (3), (4), and (7) 
> for ship roll in irregular waves — the relationships (9) and (11). 


In order to determine ship roll in irregular waves the expres- 
sion (2) and the wave energy spectrum recommended by ITTC 
was used. Results of the calculations are graphically presented 
in Fig.10 and 11. 


Explanation to Figs 10,11 : 
+ exact calculations by SEAWAY software 


approximation by using (11) 
approximation by using (2) and (7) 
= approximation by using (2, 3, 4) 
approximation by using (9) 


V [kts] 


Fig. 10. The significant ship roll amplitudes $173 , 
at V=var,d=8m,GM=0.9m,H=4m,T=10s 


From the diagrams given in Fig.10 and 11 it results that 
both presented approaches show similar interpolation capabi- 
lities. The extrapolation of roll amplitudes in irrerular waves 
by using (9) and (11) yields rather accurate results against the 


standard values. Whereas the application of the approximation 
(2) to determine significant roll amplitudes in irregular waves 
brought in erronous solutions, especially in the case of extra- 
polation. 


V [kts] 


Fig. 11. The significant ship roll amplitudes $173, 
atV=var,d=9.5m,GM=2m,H=6m,T=10s 


RECAPITULATION 


O In the paper were presented approximations of ship roll in 
regular and irregular waves by means of artificial neural 
networks and linear and non-linear regression methods, on 
the basis of main ship service parameters (speed, draught, 
initial transverse metacentric height), and — in the case of 
approximation ship roll in irregular waves — wave parame- 
ters (significant height and characteristic period of waves). 
The example approximations were elaborated for S-175 
model containership. 


O From the performed investigations it results that the rela- 
tionships obtained by means of MLP artificial neural net- 
works show the largest accuracy regarding interpolation and 
extrapolation of ship roll amplitudes both in regular and 
irregular waves. 


O Next, basing on the above mentioned approximations, two 
approaches to predicting ship roll amplitudes in irregular 
waves were analyzed : 


approximation of roll amplitudes in regular waves, and 
next calculation of their significant values in irregular 
waves 

¢ direct approximation of roll amplitudes in irregular 
waves. 


O From the performed investigations it results that the appro- 
ximation of ship roll angles by means of amplitude charac- 
teristics yields satisfactory results only in the case of inter- 
polation. Whereas the second approach provides relatively 
accurate solutions both in the case of interpolation and extra- 
polation. 


O The presented method makes it possible to approximate ship 
roll angle amplitudes calculated by using exact numerical 
methods, and they can be also applied to determine opti- 
mum ship voyage routes. The proposed approach may be 
also used to approximate real values of roll angle amplitu- 
des on the basis of data recorded onboard the ship during 
its service. 


NOMENCLATURE 


B - ship breadth 

Cpg - block coefficient 

d - ship draught 

F, - Froude number 

GM - initial transverse metacentric height 
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H, - wave height 5. Kukner A., Aydm M.: Influence of Design Parameters on 


L  - ship length between perpediculars Vertical Motions of Trawler Hull Forms in Head Seas, Marine 
T  - characteristic wave period Technology, Vol. 34, No.3, July 1997 

V - ship speed 6. Mesbashi E., Bertram V.: Empirical Design Formuale Using 
Zo - height of the ship gravity centre Artificial Neural Nets. 1st International EuroConference on 

V  - ship volumetric displacement Computer Applications and Information Technology in the 

à  - wave length Maritime Industries, COMPIT’2000. Potsdam. 2000 

b  - roll amplitude 7. Popper K.: The Logic of Scientific Discovery. Routledge, 14th 
œ  - wave frequency Printing. 1977 


8. Szelangiewicz T., Cepowski T.: Application of Artificial Neural 


Akcronyms Networks to Investigation of Ship Seakeeping Ability, Part 1. 

ICLL - International Convention on Load Lines Polish Maritime Research, Vol 8, No.3, September 2001 

MLP - Multilayer Perceptron 9. Szelangiewicz T., Cepowski T.: Application of Artificial Neural 

RAO - Response Amplitude Operator Networks to Investigation of Ship Seakeeping Ability, Part 2. 

RBF - Radial Basic Functions Polish Maritime Research, Vol 8, No.4, December 2001 

UBS - Universal Bulk Cargo Ship 10.Zhou, Z., Zhou, D., Xie, N: A Seakeeping Experiment Research on 
Flokstra Container Ship Model. Technical Report 4 of Study on 
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Task force 2 workshop 
of European Federation of Corrosion (EFC) 


Under auspices of EFC technical scientific projects #% Darowicki K., Orlikowski J., Arutunow A.: Detailed 

are carried out as the task force 2 dealing with : analysis of the passive-active transition during the pas- 

. X sive layer cracking using dynamic electrochemical im- 
Corrosion and protection of steel structures pedance spectroscopy 

against corrosion * Darowicki K, Slepski P.: Electrode impedance measu- 


rement of non-stationary systems 

% Darowicki K., Szociński M.: Impedance investigation 
of organic coatings subjected to alternating mechanical 
stress impact 


The activity is supervised by Prof. K. Darowicki, Head 
of the Department of Electrochemistry, Corrosion and Ma- 
terials Engineering, Chemical Faculty, Gdańsk University 


er Technology: #% Darowicki K., Zieliński A.: Application of electroche- 

On 20-21 June 2004 in Prague, the Department in co- mical noise technique in corrosion monitoring 
operation with the Department of Metals and Corrosion # Klenowicz Z., Darowicki K., Krakowiak S., Krakowiak 
Engineering, Czech's Institute of Chemical Technology, or- A.: Flow monitoring shows the way to diminish erosion 
ganized the workshop on : of tube outer surfaces 


# Krakowiak S., Darowicki K., Slepski P.: Pitting corro- 


Industrial heat exchanger problems — : . ; ; 
sion — practical experience and laboratory investiga- 


— Non-destructive Testing (NDT) inspection 


tions 
Results of altogether 17 projects were presented, 12 out ** Miszczyk A., Darowicki K.: Intercoat adhesion moni- 
of which were prepared by scientific workers from Gdańsk toring using impedance spectroscopy 
University of Technology. It was the following presenta- % Orlikowski J., Krakowiak S., Slepski P., Darowicki K., 
tions : Arutunow A.: Digital monitoring system of corrosion in 


electrochemical environments 
Zakowski K.: A time-frequency method for detection of 
electromagnetic field interference on metal constructions. 


* Darowicki K., Felisiak W.: Applications of ellipsometry x 
in corrosion measurements 
% Darowicki K., Kawula J.: Applications of conducting 


polymers to corrosion protection — a review Apart from the above, 4 reports were presented 
% Darowicki K., Mirakowski A.: Use of acoustic emission by scientists from Chech Republic 
LS to detect the pitting corrosion of aluminum and one from Portugal. J 
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OPERATION & ECONOMY 


Influence of running ship diesel engines on 
mixtures of fuel oil and rape oil methyl esters — 
— experimental tests 


Jacek Krzyzanowski 
Kazimierz Witkowski 
Gdynia Maritime University 


ABSTRACT 


The paper presents results of experimental tests being continuation of the research project described in the 

preliminary report [4]. The tests were carried out on a ship diesel engine supplied with marine diesel oil, 

and the oil and rape oil methyl esters mixed in different proportions. In the tests special attention was paid 

to influence of combustion of the mixtures on exhaust gas content, including its noxious components, 
as well as on possible changes of indicated working parameters of the engine. 


Keywords : tests, ship engines, renewable fuels, alternative fuels, ecology 


INTRODUCTION 


In [4] the authors presented an analysis of possible applica- 
tion of alternative fuels to ship diesel engines. As it results 
from the subject-matter literature on the state of research in 
this area, for diesel engines possible application of mixtures of 
diesel oil and vegetable oil esters is most seriously considered. 
In European conditions it first of all concerns the mixture con- 
taining from a few to a dozen or so percent of rape oil esters 
(RME). Though many tests have been performed, their results 
are often contradictory, as they have been carried out in very 
different conditions (different types of diesel engines tested 
under different loading). Therefore has been emphasized ne- 
cessity of carrying out tests under constant rotational speed. 
Simultaneously no report concerning application of alternati- 
ve fuels in question to ship diesel engines can be found. For 
this reason the authors were encouraged to undertake the labo- 
ratory tests satisfying the above mentioned premises. 


LABORATORY TESTS 
Object of tests 


The tests were carried out with the use of the one-cylinder, 
two-stroke, crosshead supercharged diesel engine which is an 
element of the test stand already used for the previous investi- 
gations on emission of exhaust gas noxious components, pre- 
sented by the authors in [1, 2, 3]. 

The test stand makes it possible to load the engine both 
with torque and rotational speed. During operation of the engi- 
ne its most important working parameters, including those elec- 
tronically indicated, can be recorded. An applied analyzer al- 
lows to investigate exhaust gas content. 


To supply the engine during the tests in question the mari- 
ne diesel oil (MDO) and its mixtures with rape oil esters (RME) 
of the following proportions, were prepared : 


> 5% of RME in MDO 
> 10% of RME in MDO. 


The MDO was of the density of 831 kg/m’, and the RME 
of 883 kg/m’. As a result of the mixing the biofuel of the den- 
sity of : 833 kg/m? in the first case, and of 836 kg/m? in the 
second case, was obtained. Detail properties of the marine die- 
sel oil and the above specified mixtures were given in [4]. 


Test program 


The tests were carried out within the broad range of engi- 
ne’s loading, namely : at 25, 40, 50, 60, 70, 80 % M/M; , and 
for constant rotational speed of the engine, set twofold : at 220 
rpm and 320 rpm. 

Ata given rotational speed and successively set loads, were 
realized measurements of the engine’s working parameters and 
its exhaust gas content during combusting by the engine : the 
MDO alone, and the two above specified mixtures (i.e. 5 % 
RME in MDO, and 10 % RME in MDO). 

The results obtained from the tests when supplying the en- 
gine with the MDO alone was assumed the reference point for 
determination of influence of combustion of the MDO/RME 
mixtures on the engine’s working parameters and its exhaust 
gas content. 


Test results and their analysis 


The test results are presented in Tab.1 and 2, and the chan- 
ges of values of selected working parameters of the engine, 
and of exhaust gas content are graphically shown in Fig.1+10 
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Tab. 1. Test results — exhaust gas content and values of selected working 
parameters of the engine in function of loading level and kind of fuel, 
at the constant engine speed of 220 rpm 
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for different loads, three kinds of the applied fuels, and a given 
rotational speed set constant. 

Analysis of the results indicates a noticable influence of 
combustion of the used MDO/RME mixtures on the engine’s 
working parameters and its exhaust gas content against those 
obtained during combusting the MDO alone. 

Within the entire range of the set loads and for both rotatio- 
nal speeds of the engine (220 and 320 rpm) a small drop of the 
maximum combustion pressure Pmax , namely by about 3% 
only (see Fig. 1 and 2), can be observed. Simultaneously, a small 
increase of the respective angle of occurrence of pmax , measu- 
red from the top dead centre (TDC) of engine’s piston, (in Tab. 1 
and 2 this is the quantity pmax), can be observed. The pheno- 
menon may reveal a somewhat longer time of combustion pro- 
cess of the applied biofuels against that of the MDO alone. 

Also, the mean indicated pressure p; drops to a small de- 
gree only. A more distinct drop, by about 3% on average, can 
be observed during running the engine on the mixture (MDO + 
+ 10%RME) at the engine’s speed of 220 rpm, and by about 
4% on average at the speed of 320 rpm (see Fig.3 and 4). The 
small drop of p; resulted also in a small increase of the specific 
fuel oil consumption fe. However in some cases the changes 
exceeded the measurement error limits only insignificantly. 

During the tests the course of pressure in the engine’s in- 
jection system was recorded. Regardless of a kind of fuel, chan- 
ges of the pressure for both the engine’s speed were small. It is 
not possible to state a.o. any significant influence of the tested 
fuels on change of the maximum injection pressure Pmax.in 
(see Fig.5 and 6). 
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MDO +10 % RME 
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Tab.2. Test results — exhaust gas content and values of selected working 
parameters of the engine in function of loading level and kind of fuel, 
at the constant engine speed of 320 rpm 

_ Values of peered 


Loading c 
at [mm | 0, [c0] no, [OO] [Pantanal C 
eec aa e 
ira] 121 [305 a |26 fozsalas7] 60 [2.6] 20 
[so fnz| ae 92 oes o0 [na 
Fao fraf se [2 ose] os | 
[ase s| 16s para 
[as [nsf] 1 [27 pz fr] ex [>| 
[a> fnsh nee] [2 fos] [v9] 
[a frf is] 10 pass] 17 >| 
[_[a|5| |> oar n 
[9a [> omens nso 
[rs || | oe] ons ane 
[sw ||| os] na 


On the basis of the exhaust gas analysis it can be stated that 
combustion of MDO with 5% addition of RME caused on ave- 
rage the drop of NO, content by over 6%, and that during com- 
busting the MDO with 10% addition of RME the drop on ave- 
rage exceeded 8% (see Fig.7 and 8). However ambiguous are 
changes of CO content in exhaust gas. As, depending on engi- 
ne’s load, both either a distinct increase or a drop of CO con- 
tent in exhaust gas, can be observed. For instance, in the case 
of using both the mixtures at 220 rpm engine speed the drop of 
CO content in exhaust gas exceeded 15% at the load levels of 
40 and 70% M/M, , but at the load levels of 25, 50, 60 and 80% 
M/M, the increase of CO content in exhaust gas by a few to 
a dozen or so percent see Fig.9 and 10) was observed. 


MDO +5 % RME 


MDO +10 % RME 


Explanation to Figs 1+10: 


m MDO AS5%RME X10% RME 
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Rys. 1. Maximum combustion pressure Pmax in function of engine load for 
different kinds of fuel at constant engine rotational speed n = 220 obr/min 
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Rys. 2. Maximum combustion pressure Pmax in function of engine load for 
different kinds of fuel at constant engine rotational speed n = 320 obr/min 
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Rys. 3. Maximum mean indicated pressure p; in function of engine load for 
different kinds of fuel at constant engine rotational speed n = 220 obr/min 


Rys. 4. Maximum mean indicated pressure p; in function of engine load for 
different kinds of fuel at constant engine rotational speed n = 320 obr/min 
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Rys. 5. Maximum fuel injection pressure Pmax, in in function of engine load for 
different kinds of fuel at constant engine rotational speed n = 220 obr/min 
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Rys. 6. Maximum fuel injection pressure Pinay, in in function of engine load for 
different kinds of fuel at constant engine rotational speed n = 320 obr/min 
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Rys. 7. NO, content in exhaust gas in function of engine load for different 
kinds of fuel at constant engine rotational speed n = 220 obr/min 
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Rys. 8. NO, content in exhaust gas in function of engine load for different 
kinds of fuel at constant engine rotational speed n = 320 obr/min 
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Rys. 9. CO content in exhaust gas in function of engine load for different 
kinds of fuel at constant engine rotational speed n = 220 obr/min 


CO content [ppm] 


Rys. 10. CO content in exhaust gas in function of engine load for different 
kinds of fuel at constant engine rotational speed n = 320 obr/min 


CONCLUSIONS 


«* Running the engine on the MDO/RME mixture made the 
pressures Pmax and p; dropping — especially when using 
MDO/10% RME fuel — at simultaneous maintaining the 
engine’s speed and torque load constant. It shows that com- 
bustion of the fuel proceeded mildly. Most probably, the 
combustion proceeded more orderly and less dynamically 
than in the case of combustion of MDO alone. It can be 
explained by the greater cetan number of RME relative to 
that of MDO. The increase of cetan number makes the igni- 
tion-lag shorter and the work of the engine ,,soft” , i.e. at 
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a more moderate increase of combustion pressure. The de- 
scribed probable course of combustion process led to only 
a small rise of specific fuel consumption. 

The observed , however small, drop of NO, content seems 
to confirm the thesis on a more moderate course of combu- 
stion process of the fuels containing esters. It may go to 
show that the maximum combustion temperature was so- 
mewhat lower. 

The ambiguous changes of CO content in exhaust gas are 
difficult to explain in the present phase of the research. It is 
known that the excess air factor of diesel engines highly 
varies along with engine’s loading. It is favourable that for 
MDO/RME mixtures CO content in exhaust gas signifi- 
cantly drops at high engine’s loads : namely at n = 220 rpm 
it starts to occur beginning from 75%, and at n = 320 rpm 
from 70% rated load. It goes to show that the loss due to 
incomplete combustion appears lower, which indirectly 
confirms the above formulated theses. This is especially 
important for ship diesel engines which usually operate 
under 80+100% rated load. 

Lack of important differences between values of the maxi- 
mum fuel injection pressures should be justifed positively, 
as it shows that a little greater viscosity of RME against 
that of MDO does not detrimentally influence operation of 
fuel injectors. 

In the light of the above presented analysis and resulting 
conclusions it seems reasonable to continue the tests with 
the MDO/RME mixtures which would have even greater 
content of the esters. 


NOMENCLATURE 


CO — 


carbon monoxide 


CO,  — carbon dioxide 

° CSR — crank angle [deg] 

fe — specific fuel consumption 
M — set torque of engine 


— 


C gines. 
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C onferen O 
SemEko 2004 i 


For already three years the SemEko scientific meetings 
have been organized by Prof. L. Piaseczny, Head of Me- 
chanic-Electric Faculty, Polish Naval University, within 
the frame of activity of Maritime Technology Unit, Sec- 
tion of Transport Means, Transport Committee, Polish 
Academy of Sciences. 


In the last year two seminars of the kind had place 
during which the following papers were presented : 


* Camless electro-magnetic timing gear for four-stroke 
combustion engine — by K. Zbierski (Lodz University 
of Technology) 

Optimization problems of control of piston combustion 
engines — by Z. Chtopek (Warsaw University of Tech- 
nology) 

and two papers prepared by J. Pielecha 
(Poznan University of Technology) : 
* Investigations of emission of noxious components con- 
tained in exhaust gas during starting the engine 

* On possible lowering NO, emission from diesel en- 


a 
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M, — rated torque of engine 

NO, — nitrogen oxides 

O2 — oxygen 

P,(MIP) — mean indicated pressure 

Pmax  — Maximum combustion pressure 

Pinaxin — Maximum fuel injection pressure 

TDC — Top Dead Centre of engine piston 

Oy max — angle of Pmax occurrence — measured as crank angle 

relative to TDC 
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Z IM Closing the 2004-year activity 


On 4 November 2004 members of the Marine Techno- 
logy Unit (acting within the Transport Technical Means 
Section, Transport Committee, Polish Academy of Scien- 
ces) met in Maritime University of Szczecin to held its 
plenary scientific session. 


During the scientific part 
of the meeting two papers were presented : 


% Ship routing on oceans — by B.Wisniewski (Maritime 
University of Szczecin) 

* = Ultimate strength analysis of ship hull — by M.Taczata 
(Technical University of Szczecin) 


After discussion on the presented papers the Unit's 
members adopted the report on the Unit's activity in 2004, 
presented by Prof. Jerzy Girtler (Gdansk University of 
Technology), the Chairman of the Unit, as well as some 
proposals to the program of the Unit's activities in 2005 
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